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Preface 


The presented material was created within the Erasmus+ project of the student Stefanos Syllignakis under the 
leadership of Petr Vosynek. It is basically support material for the subject 6KP and its English version 6KP-A (basics of 
computational modeling using finite element method) taught in the Institute of Solid Mechanics, Mechatronics and 
Biomechanics, Faculty of Mechanical Engineering, Brno University of Technology. 

‘Computer labs of 6KP and 6KP-A are composed of active exercises under the current interpretation of the 
fundamentals associated with the type of elements and also from a separate project for a group of students. The texts were 
made in two versions, for the computing open_source system Salome_Meca (C_A) and for computing system ANSYS 
Workbench v16.2. 
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Introductory 


General Principles of Code and Platform 


# General Principles of Code and Platform: 
© Code_Aster is a “stand-alone” thermos-mechanical solver. 
= No integrated GUI to create geometries and meshes. 
* No colorful post-processing 
* With study data prepared in a text file. 
= Input: Mesh and data setting. 
= Output: Physical fields (displacement, stress-strain, temperature). 


Post-process! 


Creation of finite of results 
element model and 


mechanical 


# Salome is a generic framework for pre- and post-processing. 


Various ] Job 
physical | Code Aster) Solver X Solver Y 
solvers manager 
] ] 
) ) 


Supervision, GUI and SDM 


Sal 
gener’ | CAD Data Math: 
services setting operations 
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# = Code_Aster is a thermo-mechanical solver. 


# Salome is a generic platform for physical simulations. 
© Salome_Meca = Salome + Code_Aster (Software integration). 


co 
# Advantages: 


© Easy installation of a complete framework (Linux). 
© Acconsistent and continuous graphical environment. 
= Access from different modules to main Salome study elements: meshes, results. 
= Graphical selection of topological entities for data setting of Code_Aster. 
© Possibility of using different pre- and post-processing tools. 
= Importation of meshes and geometries prepared by GEOM and SMESH Salome modules. 
= Importation of different input mesh formats and output result formats. 
© Possibility for a “stand-alone” use of Code_Aster solver. 


# Salome_Meca is the integration of Code_Aster solver in the Salome Platform. 


File Edit View Tools Window Help 


JOa@a xlee 


Geometry Module. 
Post-Pro Module. 
ParaVis Module. 


| “Salome Meca 


YACS Module. 
JobManager. 


OpenTurns. Eficas Module. 
Aster Module. Homard. 
Europlexus. MED. 


Presentation of Code_Aster 


# Code_Aster is an acronym for analysis of structured and thermos mechanics for studies and research, is a general 
Finite Element Analysis software coming from EDF (Electricite De France) R&D Department, 


#  Itcan handle mechanical thermal and associated phenome in all sort of analysis: Linear Statics, Non-Linear 
Statics, Dynamics, Thermics and more. 
# Due to its numerous capabilities, Code_Aster is a very complex affair and its somewhat unfriendly user interface 
makes the learning curve quite steep at the beginning. 
# Anall-purpose code for thermos-mechanical study of structures. 
© With a wide variety of models: 
= More than 400 finite elements: 3D, 2D, Shells, Beams, Pipes. 
= More than 100 constitutive laws. 
= A wide range of solvers: Mechanical Statics and Dynamics, Vibrations, Modal and Harmonic 
Analysis, Thermo-Hydro-Mechanical Coupled Problems, Thermic, Acoustics, etc. 
© A computational software used by engineers, experts and researchers. 
= Studies: A need of robust, reliable, tested and qualified industrial simulation code at EDF. 
= Researches: Continuous integration of new models in the development versions. 


# Solving three types of non-linear problems, 
© Material Behavior: About a hundred nonlinear constitutive laws. 
© Kinematics: Large displacements, large strains, large rotations, 
© Contact and/ or friction. 
# Advanced features in mechanics. 
© Porous media, fracture mechanics, fatigue, damage, metallurgy, seismic analysis, rotating 
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# Code_Aster Definitions. 
© FONCTIONS [FR]: Tabulated (discrete) function depending on one parameter. 
© NAPPE [FR]: Tabulated (discrete) function depending on two parameters. 
© FORMULE [FR}: Continuum formula depending on several parameters. 

# Code_Aster Parameters: 
© ABSC_CURV > Curvilinear abscissa. 

DX, DY, DZ > Displacements along X, Y, Z. 

DRX, DRY,DRZ > Rotation along X, Y, Z. 

> Coordinates X, Y, Z. 

> Strain. 

> Stress. 

> Time. 

> Temperature. 


> Command File, 

> MED file containing mesh. 
> Text file containing output of the solver (message, errors, warnings, ete.). 
> Text file containing results of the simulation in a table format. 

> MED file containing results of the simulation. 

> Folder containing the output database of the solver. 

> File containing the base of calculations in HDF format. 

> Launch script generated by the service. 


#  Eficas standard concepts: 
o DEBUT [Start]. 
Material definition, 
‘What type of Mesh to read. 
‘Type of modeling [1D, 2D, 3D, mechanical, thermal]. 
Assign previously defined material to the model. 
Add geometric boundary conditions. 
Add load boundary conditions. 
Define the type of analysis [static, linear, non-linear, dynamic, ete. ]. 
Results to be calculated at the elements and nodes. 
Save the results to a MED file. 
FIN [Finish]. 


0000000000 


Wide Range of Finite Elements 


# Continuous mechanics. 
© 3D--> Linear, quadratic, reduced or full integration, 
© 2D --> Plane strain, plane stress, axi-symmetry. 
--> Integration of non-linear behavior in plane stress. 
# Structural elements. 
© 2Delements: Shells, plates 
© IDelements: Beams, bars, cables, pipes... 


# Connections and assemblies. 

© Linear relationships between degrees of freedom, transmission of torques effort, 
# Discontinuous media (cracks and joints). 

o XFEM level-sets. 

© Joint elements and CZM (Cohesive Zone Model). 
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3k 
Algorithms and Analysis Methods 


# Mechanical solvers. 
© Linear or non-linear static 
© Dynamic on phys 
© Modal analysis: 
© Dynamic on modal 


MECA_STATIQUE, STAT_NON_LI 
DYNA_LINE_TRAN, DYNA_NON_LI 
CALC_MODAL, MODE_ITER_*... 

DYNA_TRAN_MODAL, DYNA_VIBR. 


# Other physics. 
© Thermie: THER_LINEAIRE, THER_NON_LINE... 
© Acoustic: PHENOMENE ACOUSTIQU! 
o Metallurgy (for welding applications). 
FSI: Fluid structure interaction. 
© Thermos-hydro-mechanical coupling. 


# Tools for resolution. 

© Sub structuring, control of the nonlinear algorithms. 

© Several algebraic solvers, sequential or parallel, direct or iterative. 

© Post-processing tools: CALC_CHAMP, POST_CHAMP, POST_DYNA_* 
th Outcome > Adjusting Code_Aster > C_Aster Solver > Po: 


Sketching Geometry > Arranging Mesh Inputs > Visualizing Me: 


Processing Visual, 


Code_Aster 
solving 
‘GNUFDL cence 


CAD Modeler/ Salome (GEOM) > Mesh tool/ Salome (SMESH) > Data setting/ Eficas, Wizard > Computation /ASTK > Visualization App/ PARAVIS 


Creation of a Command File 


ficter Joos Window Helo 


# With a Wizard. “| add study case 


© Available in Salome_Meca only. Gurnee & 
© 4 different wizards in Salome_Meca 2013.1. bois 
= Linear Elastic Analysis. %, Linear elastic 
* Modal Analysis. “&, Modal analysis 
* Linear Thermal Analysis. fr Linear thermic 
* Crack Analysis. “a Crack analysis IXEM) 


# With the graphical command file editor (Eficas App). 
© Provided with Salome_Meca. 
© Provided with Code_Aster standalone. 
© Cannot handle python control flow instructions (if, while, [f = aa 
for). 
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# With your favorite text editor (ex. Gedit, 


Data Setting and Command Lang 


ge 


# Command file is also a Python script. 
© However, we should only focus on the Code_Aster 
commands. 
# Command file composes of a sequence of Code_Aster 
specific commands. 
© Each specific command composes of keywords and 
defines, assigns or uses data as input. 


# Most of the commands produce “concepts”. 
© On the left side of the equal sign ( 
© The concepts generated by one command can be used 

as an input to the following command, 


# Command file contains no geometry description. 


# Examples: 


= DEFI_MATERIAU (ELAS=_F(E=210000, NU=0.3, RHO=0.007850,) .); 


ra Command Pa Pi of the Simple Keyword 
Factor Keyword 


Name of the concept, User choice Simple Keyword 


co IMPR_RESU (RESU=_F (RESULTAT = RESU), ): 


re Factor Keyword \ i oe 


Command Argument of the Simple Keyword 
No Concept is Produced Simple Keyword 
° ” keyword is used to extend an existing concept. 


MATLLAGE = mes 
CREA_GROUP_NO <_}) (GROUP_MA ="upper’, ),): 


Reuse of the MA Concept Defined by LIRE_MAILLAGE oe 
Indicator of a Factor Keyword 


© Since commands are sequential, a concept must be created before being used. 
MESH = LIRE_MAILLAGE (); 
STEEL = DEFI_MATERIAU (ELAS = _F (E=20500e06, NU=0.3,) .); 
CHAMT = AFFE_MATERIAU (MAILLAGE=MESH, 

our, 

yi 
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Salome’s Possibilities 


# What can we do with Salome: 
© Import and export, repair and clean, create and modify geometries (CAD). 

Mesh, quality control, import/ export meshes. 

Manipulate the physical and numerical properties of a geometry. 

Manage the various stages of a computation: receive data, setting up a calculation, return results. 

Run the computation sequences and coupling between solvers. 

View and post-process the results. 


# Salome, GEOM Module: 


oo0000 


© Design of geometric obje 
© Import/ export of objects of different CAD formats. 
© Repairs and correction of CAD models. 

© Adaptation of CAD models for computer simulation. 
© Based on the Open Cascade technology. 


# Salome, GEOM Module: Import & Viewers 
© Various import/ export formats: STEP, IGES, BREP, STL, ACIS. 
© Graphic functionalities: Transparency, coloring, shading/ wireframe, rotation, zoom, pan. 


# Salome, GEOM Module: GUI 


File Edit View New Entity Operations Repair Measures Tbols Window Help 
Dee X|B eee le 2 Bm @ % mm em Gl 1 
sim~aeotuc|:joaf@nnd:|vosanzzasey? 
Awe saalsaoagePieauviavelen [was a va « 


7 
‘4 2eCOCl: Bae ain ay ale + ow © 
2 
@ > Salome Components, @ > Basic Primitives. @ > 3D Primitives. 
(@ > Advanced Primitives. ©) > Boolean Operations. © > Generations. 
@ > Transformations. @® > Constructions. @ > Measures. 
# Salome, GEOM Module: Terminology 

o. Vertex. 

o Edge. (Warten 

© Wire (set of edges). —, 

© Face (constructed from a closed . 


wire). 
Shell (set of 
Solid (constructed with a second 
shell). 
© Compound (set of any types). 
# Salome, GEOM Module: Conception 
© Geometric Primitiv 
= ID: Line, circle, ellipse, are, curve, vector. 
= 2D: Plan, working plan. 
= 3D: Parallelepiped, cylinder, sphere, torus, 
* Sketch 2D/3D: Construction of complex lines or surfaces. 
* Construction of elementary objects: Vertex, edge, wire, shell, solid, compound. 
= Advanced Primitives: Tee pipe. 


cone. 
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° 


° 


° 


Operations 

= Boolean: Merge, join, cut, intersect. 

* Fillet, chamfer, partition. 

= Explode: decomposition into basic objects. 
Generations 

* Extrusion, revolution, filling, pipes. 
‘Transformations 

+ Translate, rotate, mirror, scale, offset. 

= Multi-translation, multi-rotation, 


# Salome, GEOM Module: Partitions and Objects 


° 


° 


° 


Define the relevant topological entities to ease mesh creation, calculation setup and post-processing 
result. 

Geometric Groups: Manage geometric objects from which can be created element/ node groups required 
when assigning boundary conditions or material properties. 

Geometric Operations to Define Partitions in a Mesh: Non-manifold geometry, for hexahedral meshing 
splitting non-hexahedral shape into hexahedral ones, define different refinements, merge. 


# Salome, MESH Module 


° 


° 


The geometric model is imported from GEOM Module. 

Information and quality control of meshes. 

Groups of nodes or elements and operations on these groups. 

Various import/export formats: MED, UNV, STL, CHNS. 

Modifications of meshes 

Principles: 
* One algorithm for each dimension [1D: Wire Discretization]. 
* One hypothesis for each algorithm [Nb. Segments]. 

The mesh calculation starts from the smallest to the largest dimension: 
* 0D91D>2D>3D. 


# Salome, SMESH Module: Algorithms 


° 


° 


For each dimension, several algorithms are available: 
= 1D: Wire Discretization, Projection 1D... 
= 2D: Quadrangle (mapping), NETGEN 1D-2D, BLSurf, 
* 3D: Tetrahedrons (NETGEN), Tetrahedrons (GHS3D), Hexahedrons (i, j, k)... 
Some algorithms are multidimensional: 
= 1D-2D: BLSurf, NETGEN ID-2D... 
*  1D-2D-3D: NETGEN ID-2D-3D... 


# Salome, SMESH Module: Hypothesis 


° 
° 


° 


‘The parameters of a meshing algorithm are set through a Hypothesis. 
Ex, for 1D mesh: 

* Algorithm 1D: Wire Discretization. 

= Associated Hypothesis: Nb. Segments = 4, Equidistant Distribution. 
Ex, for 2D mesh: 

= Algorithm 1D: Wire Discretization. 

= Associated Hypothesis: Nb. Segments = 4, Equidistant Distribution 

* Algorithm 2D: Quadrangle (mapping). 


# Salome, PARAVIS Module: Terminology 


° 


The Salome post-processing module based on ParaView. 
= With functionalities added by EDF. 
© MED interface. 
© Integration points. 
* Modal animation, 


(ag Ntrobuction FINITE ELEMENT ANALYSIS METHOD 


o Filters. 
= In ParaViS, the data is managed by means of filters. 
* Ex. Deformed Shape, cutting plane ... 
= The filters depend on the type of data. They can be chained. 
© Initial data > Filters 1 > Filter 2 >... 


o Views. 
* A surface on which you can see the data. 
* Ex. 3D, 2D, histogram, plot... 
© Displays. 
=” The data can be seen in different ways in different views. 
* Ex, Surface, Wireframe, PointSprite... 


Code_Aster Possibilities 


# Special commands: DEBUT, FIN, POURSUITE, 
© The DEBUT command: 
= Begins execution, ignoring the previous lines. 
= Sets the logical units to be used for messages and files results. 
= Basic usage: DEBUT () 
© The POURSUITE command: 
= Restarts execution from a base provided as an input. 
* Useful to continue a calculation initiated with the same version of Code_Aster. 
= Recommended to decouple the calculation of post-treatment. 
= Basic usage: Calculation: =DEBUT() —.... FIN() 


Post: POURSUITE()...  FIN() 


© The FIN command: 
= Ends the command file and ends the run, ignoring the following lines. 
= Closes the base at the end of execution: folder containing all the concepts generated during the 
calculation (mesh data, intermediate structures, results). 
* Specifies the format used for the produced base: .hdf or Aster, 


# What is a Mesh? 
© Coordinates of nodes. 
© Cells defined by their connectivity. 
© Groups of cells (GROUP_MA) and groups of nodes (GROUP_NO). 


eater eS ra 


Nodes Cells GROUP_MA GROUP_NO 
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# = Read the Mesh 
© Meshes at Aster format or MED formats: 
= Code_Aster format: 
* mymesh = LIRE_MAILLAGE () 
* MED format: 
© mymesh = LIRE_MAILLAGE (FORMA’ 
© Meshes of other formats: 
* Commands PRE_***: PRE_GIBI, PRE_IDEAS, PRE_GMSH. 
= Ex. PRE_**() 


‘MED’ ) 


mymesh = LIRE_MAILLAGE () 


o Results: 
* Command IMPR_RESU, 
* For instance in MED format: IMPR_RESU ( FORMAT = ‘MED’, 


RESU =_F (RESULTAT = myresult,) ,); 


# The Choice of Finite Elements 
oA finite element is: 
* A geometric description, provided by the mesh. 
+ Shape functions. 
* Degrees of freedom. 
© Its choice determines: 
* The equations that are solved. 
+ Discretization and integration hypothesis. 
* The output fields. 
© 3D elements examples: 
* MODELISATIO! 


© 2D element examples (plane strain, plane stress, axi-symmetry): 
* MODELISATION = ‘D_PLAN? 
‘C_PLAN? 
‘AXIS’ 
© 1D element examples (beams, pipes, cables): 
* MODELISATION = ‘POU_D_T” 
‘TUYAU" 
© Shell elements examples (plane strain, plane stress, axi-symmetry): 
* MODELISATION= ‘DKT? 
‘DST’ 
‘COQUE_3D° 
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# The Definition and Assignment of Materials: 
oA material is the definition of numerical parameters 


Units of measurement are not managed by Code_Aster, the user must use a consistent system of 
units between materials, unit length of the mesh, loading data, like: 
* Coordinates of the mesh nodes: mm > m. 
© Young’s modulus: MPa > Pa. 
Applied force: N > N. 
© Stress obtained as a result: MPa > Pa. 


© The material must be assigned. 


Assignment to a geometrical area: Groups of finite elements designated by the name of groups of 
cells. 


© The material must be consistent with the assumptions of resolution (constitutive law). 
© Be careful to assign what has been defined, then use what has been assigned. 


# = The Characteris 


's of Structural Elements: 


© Shells, plates, beam, pipes, discrete elements... 
© Geometric information not provided by the mesh: 


# Loadings: 


Shells: Thickness, reference direction in the tangent plane. 

Beams: Cross-section, orientation of the principal axes of inertia around the neutral axis, 
curvature of the curved elements, 

Discrete Elements (spring, mass/ inertia, damping): Values of the stiffness, mass or damping 
matrices. 

Bars or Cables: Area of Cross-section, 

3D or 2D Continuous Media Elements: Local axes defined with respect to the anisotropy 
directions. 


© Commands AFFE_CHAR_*** 
© Loadings inside the continuous media: PESANTEUR, ROTATION, FORCE_NODALE... 
o Edge loadings: FORCE_FACE, FORCE_ARETE, PRES_REP... 


oo 


Loadings specific to structural elements: FORCE_POUTRE, FORCE_COQUE, FORCE_TUY. 
Pre-established relationships on the degrees of freedom: 


DDL_IMPO: DoF imposed on a node or group. 
FACE_IMPO: DoF imposed on the nodes of a cell or a group of cells. 
LIAISON_SOLIDE: Regasification of a set of nodes or cells. 


© Be careful to the constistency with the degrees of freedom allowed by the chosen finite element. 


# Solving the Problem: 
© There are about fifteen resolution operators to solve the physical problems: 


Thermic: THER_LINEAIRE, THER_NON_LINE 
Mechanics: MECA_STATIQUE, STAT_NON_LINE 

Dynamics: DYNA_VIBRA, DYNA_NON_LINE 

Modal Calculation: MODE_ITER_SIMULT, MODE_ITER_INV 


© One must give the description of the problem set prior solving. 


Materials: CHAM_MATER 
Geometrical Characteristics: CARA_ELEM 
Loading: EXCIT 


A time list if necessary: LIST_INST 


© One can also change settings on the resolution algorithm, 


In most cases, the default values are suitable. 
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# Choosing the Linear Solver. 
©. The solver can be chosen via the keyword SOLVEUR/ METHODE. 
© Relevant depending on the problem to be solved. 
o Direct solvers: 
* MULT_FRONT (multi-frontal): Default method. Universal solver, yery robust. Not 
recommended for mixed models of X-FEM, incompressible... 
= MUMPS (multi-frontal massively parallel sparse direct solver): External solver. Same scope as 
MULT_FRONT. Provides access to parallelism. 
© Hterative solvers: 
+ GCPC (preconditioned conjugate gradient): Recommended method for thermic. Useful for well- 
conditioned, large problems. 
+ PETSC: External multi-method solver. Very robust when associated with pre-conditioner. 


# Fields of Quantity. 
© Location of the values: 
"Fields on nodes: NOEU; 
= On the elements: 
* Field by element on Gauss Points: ELG 
* Field by element on Node 
© Constant field on element: E 
o Naming rules XXXX_YYYY: 
= 4 first letters: name of the quantity, SIGM, EPSI, ERRE; 
* 4 middle letters: location, NOEU, ELGA, ELNO, ELEM; 
© Location of the values: 
* Three locations: 


_ELNO NoEU 


_ ~ 
CHAM_ELEM CHAM _ELEM cCHAM_NO 
given at Gauss points given at nodes. given atnodes 
SIEF_ELGA per cell 
a srr TNO. STEF_NOEU 


# Concept Resultat. 
© Resolution operators produce typed data structures: resultat. 
© The type depends on the operator: EVOL_ELAS, EVOL_NOLI, EVOL_THER, MODE_MECA.... 
© Ateach computation step, one or more fields are stored in the data structure: resultat. 
© Fields are identified by access variables: INST, NUME_ORDRE, FREQ, NUME_MODE... 
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# Calculate Physical Fields or Perform Mathematical Operation: CALC_CHAMP /POST_ELEM. 
© Produce fields and tables. 
o CALC_CHAMP 
= Create or complete a DS resultat. 
= Calculation of stresses, strains, energies, criteria, error indicators... 
o POST_ELEM 
= Create a table. 
* Calculation of energies (potential, elastic, kinetic, dissipated..). 
= Calculation of integrals or averages. 
# Handling Fields and Tables: CREA_CHAMP/ CREA_RESU/ CREA_TABLE/ CALC_TABLE. 
o CREA_CHAMP: Create fields of quantitates champ_gd. 
o CREA_RESU: Create a DS resultat from fields, 
o CREA_TABLE: Create a table from a function or a list of real. 
o CALC_TABLE: Manipulate data’s from tables like a spreadsheet. 
# Print a Result or Fields: IMPR_RESU. 
o Print: 
= Meshes. 
= Fields of quantities. 
= Data contained in a DS. 


© Different formats: 
= Listing: *‘RESULTAT” and ‘ASTER’ 
"| MED: ‘MED’ 
" GMSH MSH" 
= CASTEM *CASTEM’ 
= IDEAS FORMAT = ‘IDEAS’ 


# = Print a Table. 

o Print the contents of a table in a listing or an Excel file. 

© Also plot curves. 

o Selection: 
* NOM_PARA: Choice of columns to print. 
+ FILTRE: Choice of lines to print (some lines checking criteria). 
+ ‘PRI: Choice of the order of rows to print (ascending or descending). 
* FORMAT: Choice of print format. 


Where to Find Information 


© U documents: Usage. 
= Using the commands. 
= Methodological documents, 
* Documents for starting: 
© [U1.00,00] Getting started with Code_Aster 
© [U2.02.00] Introduction to Code_Aster 
© V documents: Validation, 
= One test-case = one V document. 
* Description of the problem and the reference solution. 
© Rdocuments: Reference. 
= Theory related to methods and algorithms of the code. 
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# Tools for the community. 
© The forum is: 
= A place of exchange and mutual assistance for Code_Aster users. 
= Exchange of best practices, tips. Sharing of experience. 


Useful Links 
# www.code-aster.org > Official Code_Aster WebSite. 
#  www.salome-platform.org ~ Official Salome Platform WebSite. 
#  hup://code-aster.org/forum2/ > Official Code_Aster, Salome Meca Forums. 
# — http://www.code-aster.org/spip.php?rubrique19 > Code_Aster, Salome_Meca Documentation. 
# http://www.code-aster.org/spip.php?rubrique68_ > Code_Aster, Salome_Meca Training. 
# http://www.code-aster.org/spip.php?rubrique22_ -> Code_Aster, Salome_Meca Support. 
# — http://caelinux.com/CMS/index.php?option=com_content&-view=article&id=55:caelinux- 
2013é&catid=694ltemid=41 > Download CAELinux, Salome_Meca WebSite. 
#  https:/Avww.virtualbox.org/ > VirtualBox Download WebSite. 


# —hup://caelinux.org/wiki/index.php/Doc:CAETutorials > Code_Aster, Salome_Meca Tutorials, 


Personal Difficulties During Salome_Meca, Code_Aster Study 


# The first problem I faced concerns Salome_Meca Post-Processing Module. I noticed that even though my mesh 
and my command file where absolutely correct, the platform was running the analysis BUT Post-Pro branch 
didn’t appear after the analysis. I figured later on, that this was a problem concerning the names I was giving to 
the folder, files, command files, aster study. That’s why you should avoid giving difficult names or even adding a 
space between letters, stay simple and also avoid symbols. 


# Secondly, Geometry Viewer mode will in most cases get stuck. You will need to reset it everytime that happens, I 
found a way by closing the Geometry viewer, activating Mesh Module and then re-activate Geometry Module 
without making any changes. 

# Lines beginning with “#” are comments. Also, many things are explained directly in the code with commented 
lines at the right place. A line with [U4.81,03] for example, means that the following command or operator is 
described in the U4.81.03 document in Code_Aster documentation available at the links given above. 


# Another important note, is that within the command file we must not use concept names more than 8 characters. 
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Chapter_I: Cube 


In the first chapter we will get familiar with the software concepts 
and in general with the idea of Finite Element Analysis Method 
We will work on a Cub 


-ometry to easy our prelusion 


Step 1: Purpose of the FE Analysis/ Description of the problem 


The purpos 
getting to know the 


of this Finite Element Analysis is to make an outset of the Fini 
ncepts of our Open_Source software, Salome_Meca. 


Element Analysis Method, and also 


Step 2: Input Values for the FE Analysis 


Cube Dimension; 100mm x 100mm x 100mm; 
Young's Modulus: 2.1e5 MPa; 
Poisson's Ratio: 0.35 


sure Applied on the “Top” face: i) 5 MP: 
ii) 50 MPa; 
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Step 3: Model Geometry 


# | Description Figure 
7. | After installing Salome_Meca in your Ubuntu system we are ready to run the application. There are two different 
ways to start Salome_Mec: 
First one is to open up a terminal and go to the folder where Salome_Meca is installed and type ./runAppi 
Or double-click on Salome_Meca icon on desktop. 
When Salome_Meca loads, it will look like the figure below. 
2 TE CEI 
(Da TeU a Peae sere oueeens sane s 
3. | The figure below shows the shape of Salome_Meca toolbar. 
The first line has a pretty common menu bar: File, Edit, View, Tools, Window, Help. 
The second line, except the New, Open, Save, Close, Copy and Paste tools, consists of (1) Geometry, (2) Mesh, (3) 
Post-Pro, (4) YACS, (5) ParaViS, (6) Jobmanager, (7) Openturns, (8) Eficas, (9) Aster, ) Homard, (\)) Europleax, 
(1) MED, and a drop-down list which contains all of the above. 
rar Sle mecs 20151 
File Edit_View Tools window Help 
[Dad lad 8D & |Frscerenes lee ee FG mw me 
5. | Let’s start drawing our model. 7 ‘Activate module x 
Select Geometry either from the drop-down list 
or by clicking on the icon, STs SV eCOg WUE weer ey) 
Select “New” from the pop-up window that & Please, salect required action by pressing the corresponding button below. 
appears. 
open... | toad... | Load script. cance, 
See figure on the right 


yy Carrer: CuBe 


GeTTING_TO_KNOW_SALOME_MECA 


#_| Description 


6. | Click on New Entity > Primitives > Box from the toolbar. 


Fill in the gaps like the figure on the right. 


Finish with “Apply and Close” 


= Box Construction x 


ce ca 


(-Rosult name 


Name [0x 1 


[Aonyend ce] _aoov_| axe | wee | 


name of each button, 


even more. 


7. | Now in the viewer tab we got a new menu bar, if you move your mouse to each one of these icons it will show the 


Most of those options are used to change your geometry’s view, either by zooming in and out, rotating the model or 


(OCC scene:3 - viewer 


i 
TOW eee OE OOOO Ce Gio: VIE URY Sey 


9. | Easy as it may seem, you have 
created your geometry and it is 
time to proceed to the analysis. 


First we need to create groups 
which will be a huge help to us 
as we continue through the FE 
Analysis. 
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# | Description [Figure 1 

10. | Select Box_I from the Object Browser. > Create Group x 

shape Type- 
Got to New Entity > Group > Create from the toolbar || ae a 
and fill the gaps as shown in the figure. 
Group Name 
Choose the 3" option to select faces. ome) re 
Main Shape and Sub-shapes 

Name the group with a proper name (we used Fix 

eee proper nae (we ee mainshe @[fexa STS 
because we are targeting the fixed support face). 

Main Shape Selection estrictin shay nly selected 
Click on the bottom face to choose it and then press © No restriction ay 
“Add”. 7 Geometrical parts ofthe Second Shape | ————— 
© Only Sub-shapes ofthe Second Shane | Show al sub-shapes} 
Once you select it, it will get highlighted with white eon sens) (| | l 
edges. 
Select al 
Finish with “Apply”. ah 
Remove 

Do the same procedure but now for the face that we 
want to apply the load (top face), you can name it Load. ne =n IDES 

11. | After successfully completing all the previous steps, you should ‘Object Browser ax 
have the same Object Browser as shown on the figure on the right. 
The “eye” icons allow you to choose which part of our geometry 
will be shown on the Viewer tab. They are also used in every other 
Module (mesh, post-pro, ... ) 

12. | Always save your work before = SaveFile =| 
moving to the next step/ module. [fs romelnetmalDeakion (Course Boal Salome Mecano deooame 
Remember if you use meaningful | | Jp chnpute 
names it is going to be really iacue, 
helpful later on. fa 
Note: Be extra careful when you | | Fiename: [introductor| Save 
name your work. Do not include | | gs e¢type: 7-07 a 

space in your file names. —— 
‘Quick path: [/op1/SALOME-MECA-2013.1/SALOME-WECA-2013.1-GPLImodules/SAMPLES V6_60 >] |Add path 
# Taking in consideration that you have completed all the previous steps correctly, we will continue to the Mesh part, 
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Step 4: Meshing Geometry 
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# | Description [Figure 
7. | From the drop down list choose Mesh module. = Create mesh +x 
Name 
Select Box_1, g0 0 the menu bar Mesh > Create | | cay [famed 
Mesh. = es 
Keep the default name, Mesh_1. 30 | 20| 10 | op| 
Algorithm aa | 
Make sure Box_1 is on the geometry option and SS ——— sno 
choose Netgen 1D-2D-3D Algorithm. plete = al 
Add. Hypothesis [Hoe al 
This algorithm option will automatically fill the 3D, 
2D, 1D, 0D tabs. Asian a set of hypotheses A 
Apply and Ciose| apply | ose | Help 
2. | We need to specify a hypothesis also. ~ Hypothesie Construction +x 
5 Netgen 30 
A hypothesis isa parameter that modulates the Algorithm. In this ||" jegumante | tacts: | 
case we will adjust the Max. and Min. Size of the elements used 
by the NETGEN ID-2D-3D Algorithm. Name [RETCEN 30 Porameters 
Wax. Size | | 
Click on the gear and select Netgen 3D parameters, . Min. Size 
In the pop-up window leave everything as default (you will be Rese 
explained in next chapters about those options). 
Fineness [Moderte >] 
Click “Ok” to get back to the previous pop-up window. Growth Rate | 
Nb. Segs per Edge 15 | 
Nb. Segs perRags [I] 
F Allow Quadrangles 
F optimize 
Finish with “Apply and Close”. 
3. | Now your Object Browser should Took like the figure on the right. 


If you noticed, Mesh_1 icon has an exclamation mark, which means 
that we may have set the mesh options but we still need to compute 
it, 


To do so, right-click on the Mesh_1 and select Compute. 
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Description 


If everything is well 
placed, after some 
seconds you will get 
a pop-up window 
which contains all of 
the Mesh 
information. Like 
number of Nodes, 
Edges, Faces, 


In the figure on the 
right you can also see 
how our model will 
look after we 
compute the mesh. 


# That w 
next s 
Wizard. 


s the end of our meshing step. Remember to always save your work before going to the next step. In the 
{eps we will see how to build information for a Linear Static FE Analysis 


with the help of Salome_Meca 
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Step 5,6,7,8 and 9: Salome_Meca Linear Static Analysis Wizard 


# After completing our geometry, arranging the correct groups, using the most suitable algorithm and hypothe: 
the mesh, now is time to make some more modification for out linear static anal 


# |Description Figure 
Z. | Select Aster Module from the drop down list or by] 7 
clicking on the icon from the toolbar. | 


: ‘Qe-subapplication = 


Isotropic linear elastic study 
Mode! defintion 


Go to Aster > Wizards > Linear Elastic. 
‘what kind of mode! do you want to work on? 


Step 5: Define the type of model, select 3D like in aaa ~) 


the figure on the right. 


Proceed with “Next” 
2. | Step 6: Select the mesh that will be used for this 


specific analysis. teotropic tinea elastic study 
Tesh election 


‘Qe-subapplication x 


If for some reason Mesh_1 is not automatically } 
chosen, press the blue arrow button and then select ‘Select = mesh from the Salomé object browser 
Mesh_1 from the Object Browser. ira} 


Mesh 


“Use geometrical groups”, the ones that we defined 


© Use mesh oroups 
in Geometry Module. 


© Use geometrical groups 


Proceed with “Next” paar 
3. | Step 7: Define the material properties that we will x x 
use in our analysis. 
Wataval propertos 
Add Young’s modulus and Poisson’s ratio as 
shown in the figure on the right. Young's modulus and Poisson ratio defintions 
‘Young's modulus (E) [2.165] (E>=0) 

Our dimensions are in millimeters so our Young’s Loot) el 008) 


modulus should be in MPa. 


Proceed with “Next” < Back 


cancel 
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# |Description [Figure 
4. | Step 8: Input Geometric Boundary Conditions. 


Qt-subapplication x 


Isotropic linear elastic study 
‘This is where the geometry groups we created Boundaries conditions 


before will come in handy. 


3 ; ‘Adding imposed degrees of freedom on groups 
We want fixed support in the Fix group. fl 


‘Group ox by om 


Fixed support has zero degrees of freedom for the 
selected group (in case of volumetric geometry- 
this leads to discretization by solid elements), so 
fill in the gaps according to the figure on the right. 


Fix 0 |o ° 


Beck | Next> | cancel 
Proceed with “Next” ea | cee || Reza 


5. | Step 9: Input Load Boundary Condition, x Qesubapplicaton x | 
Isotropic finear elastic study 
Select Load and input 5 MPa if it is not already Boundaries conditions 
selected. = 


@ | asaing pressure on meshes groups 
Load is the top face group we created in the [3 


Geometry Module. ‘Group Frese 


jaftoas 


<bce cance 


Proceed with “Next” 


6. | We have defined all the information - command file selection | 
needed for the FE Analysis. Look in: [fj homejnetme/Desktap/C_ooxisalome Mecaintm =] GO © O Bi] 
neime 
We also have created our command file IS computer 
(comm), that passes on the information we |} |e cya. 
just input in the wizard, to Code_Aster. fit tn 
Click on the icon to save the (.comm) file 
in the same folder that the .hdf file is and, 
for your comfort, name it with a 
meaningful name. 
Click “Save” |) ee pame  [inmat Save 
End with “Finish” | | Flee oftve: [Cammand ile (comm) =] _ cancel 
7, | A new branch has been created to the Object Browser. 
Open “Aster” and you will see the “linear-static” underneath, 


# Summarizing our work so far, we have created our geometry, separated the boundary groups, computed the mesh and 
we created the command file containing all the necessary information about our analysis. Only thing left to do is ... 
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Step 10: Running the Analysis 


# Salome_Meca makes it easy for beginners to create and run an analysis (always depending on the kind of the 
analysis). The end results (output) can be changed based on the results we are after, by editing the command file. 
But for this introductory analysis we will use defaults, 


# | Description Figure 
T. | Right-click on “linear-static™ in the Object Browser and select "Run™ 


A new window will pop-up displaying lots of different things. 


The graphic environment is the software's way to show procedures, steps 
of solution by text output. 


The software now calculates all the information we entered before. 


2. | Ifeverything goes well, a new branch named “Post-Pro” will 
appear in the Object Browser. 


If Post-Pro branch doesn’t appear, it means that something went 6 Geometry 
wrong with the analysi > Mesh 

&- & Hypotheses 
When that happens, the information about the error is written in & & Algorithms 


text file with suffix .mess in the directory, where the task was & & Mesh_1 
saved. If you are not sure about the mistake, it would be better to 
re-do all the steps from scratch, 


There is another way of correcting those errors but we will see it in 
the next chapters. 


# Post-Pro branch is now on the Object Browser which means that the results are ready to be evaluated. Depending 
on the type of the analysis the calculation could take from a few minutes to several hours. Next step is to visualize 
the results. 
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Step 11: Post Processing the results 


GETTING_TO_KNOW_SALOME_MECA 


and select Deformed Shape. 


Open this branch, right click on “0, 


Check the Magnitude coloring in the pop-up window and keep 
the default scale factor. 


You can go through the other 2 tabs Scalar Bar, Input, but right 
now we won't need to change something else in those tabs. 


Finish with “ 


# [Description Figure 
T, | Select Post-Pro from the drop down list to enter Post-Pro Module, | Object Browser ax 
Name 
In the Object Browser we can now see three branches in the parent. |] S rcs” 
branch of field. & & Hypotheses 
= Algorithms 
ee Mesh a 
These branches contain the results of the analysis, and the only 
thing left to do is explain to the software the way to visualize those 
results, 
Groups 
S Heiss 
@) RESU___DEPL. 
~ anivades 
‘ontvades 
©. RESU__SIGM_NOEU, - 
~anNedes 
o 
2, | RESU__DEPL > Result of the Displacement of the Nodes. . Deformed shape +e 


Detames shape | seatr ase | input | 


Sealerscter: [350575 


aga rg 


‘Apply and Cose| 
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# | Description ure 
3. | A new branch appears 


inside RESU__DEPL. | © w9ynewuiroeesved cles 


aave><vonss 
With right-clicking on 
the Def Shape you get 
plenty of options. 


First of all select the 
“Show Only” option and 
then right-click > 
Representation > 
SurfaceFrame. 


‘As you can see from the 
figure at the bottom of 
the window, a graph of 
Deflection with 
minimum and maximum 
values is shown, 


4 | Now in this step, T will show you how to add more than one : Deformed Shape + 
different types of visualization in the study if required. EES] 
Scab range; 
Right-click again on “0,-“ select Deformed Shape, check the ce vals IRE Lovo 
Magnitude Box with default Scale Factor. Aiea Peas ase 
‘ Min DUO Max: [000357407 
Click on Scalar Bar tab, go to Scalar Mode and select DX from 
the given option, as shown in the figure. Colors and ab 
No, of coors: [5 = | 
‘Orientation 
© Vertical © Horizontal 
$$¢_¢_—<$<$<—<_—_——— 
» [or | 
Dimensions (n ofthe sie of View! 
wit: [5 Height: [ooe =] 
“ext properties. Bar properties. 
Preview THe salar bar 
TT show dizributon * save.as defaut 
‘esiy and Cosel [apply] _cancel_| __He 
Finish with “Apply and Close” 


CHAPTER_I: CUBE 


GETTING_TO_KNOW_SALOME_MECA 


Description 


Follow the same steps 
as before to make it 
visible and with the use 
of SurfaceFrame, 


Here you are able to see 
how the cube deforms 
in X direction. 


You can also see in Y,Z 
direction by changing 
the DX option from the 
drop down list given in 
the Scalar Bar tab. 


ure 


We can see the stress results now, and the points where 
those stresses reach maximum and minimum values. 


RESU___SIEQ_NOEU is the results (RESU) of the 
stress (SI..) Equivalent (..EQ) or equivalent stress at 
nodes (NOEU). 


Again right-click on “0,-* choose Scalar Map and from 
the pop-up window select VMIS (Von Mises), rest keep 
default 


Finish with “Ok” 


Se eC: a 


scalar Bor | input | 


Scalar range 
Sealer Node  Loganthmie 
Use fed range © Use imposed range 
Min [Eas ma: [Sasra6 
Colors an nbels 
Nb. of colors [oe =] no, of eis: [> 
Paaanaet 
© Vertical & Horizontal 
Pa 
«for —_s]¥ oor 
Dimensions (in % ofthe se of view) 
wath: [OE sigh: foce 
ee | ee 
Proview TF Hide seta bar 
F snow cistron - Save asdefaut 
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# | Description 

7. | A new branch has appeared 
now at the Object Browser, 
right-click > Show Only, 
right-click > Representation 
> SurfaceFrame to get an 
understandable view. 

Notice that at the bottom of 
the window a graph with 
minimum and maximum Von 
Mises values is shown. 

As in the deflection branch we 
can also change the fashion of 
showing (to other similar 
stresses) like TRESCA or 
YMIS_SG. All you have to do 
is follow the same steps as 
before and from the pp-up 
window instead of choosing 
MIS, select any other of the 
given options. 


Hui rueelsvyulslevuelae 


8. | Following the same steps as 


before for the “ eau 7oea)ivvuleeuseiae 
RESU___SIGM_NOEU 
branch, 


RESU__SIGM_NOEU is the 
result of the stress or SIGMa 
at Nodes. 


Here we have 6 different 
stresses, Sx, Sy, Sz, Sxy, Sxz, 
Syz. 


# Weare fi 
scenario, 


hed with our first load case, now let’s see what changes we have to make to solve our second case 
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# | Description Figure 
1. | We won't have to start from scratch since we are = Qt-subapplication x 


going to use the same geometry and mesh. sci Fa salle wap 
‘Boundaries conditions 


Choose Aster Module from the drop down window. 


Make sure Mesh_1 is selected and follow the exact |e secs ores oo we aces 


same steps we did for the first load case. 


Group, Pressure 


BUT, as it is shown in the figure on the right, when it [fons 2 
comes to input the value of the load you will have to 
input 50 MPa. 


Continue with saving with different name and finish. 


2. | As we can see from the figure, a new analysis is added | Object Browser ax 
to Aster branch. 


Right-click > Run and if everything is correct there 
will also be a new branch in Post-Pro branch. 


= Hypotheses 
Repeat the same steps as in the previous analysis to & & Algorithms 
visualize the results in this load case. & & Mesh_1 
eM Aster 


&- 3 linear-static 

= 
© @ Post-Pro 

®-* linear-static.rmed 

&* linear-staticO.rmed 


# This task was a complete FE Analysis which is very helpful to understand the steps we follow in Finite Element 
Analysis. 
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Chapter_II: Cylinder Rod 


In this next chapter we are going to get familiar with Code_Aster and how 
programming affects Salome_Meca software. 
‘We will study a Cylinder Rod case and use command line to set the necessary 
parameters, 


Step 1: Purpose of the FE Analysis/ Description of the problem 


The purpose of this Finite Element Analysis is to get one step close to understanding Finite Element Analysis 
Method, and also getting to know the concepts of our Open_Source software, Salome_Meca and Code_Aster. 


Step 2: Input Values for the FE Analysis 


Cylinder Dimensions: 1= 1000mm, r= 15mm; 
Young’s Modulus: 2e5 MPa; 
Poisson’s Ratio: 0.3; 


Pressure Applied on the “right edge” face: S00 N; 
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Step 3: Model Geometry 


EDDITING_COMMAND_FILE 


# | Description Figure 
T, | Follow the same procedure you did in the previous chapter. 
Starting from Step 3: Model Geometry #1 to #6. 
2. | Click on New Entity > Basic > Circle in Salome_Meca . ‘lrcle Construction x 
toolbar. rele 
© c e 
Fill in the gaps according to the figure on the right to create our Es S g 
Cylinder. Result name 
meme[Se 
Arguments 
Center Point (Origin by defauk)_@ | [O 
Vector [Z avs by detout] @ | loz 
Radius | 


Finish with “Apply and Close” 


Apply and Close| Apply Chose Help 


according to the figure on the right. 


Extrusion command. 


3. | Click again on New Entity > Build > Face and fill the blanks 


Finish with “Apply and Close” 


*We made a face from our geometry, to be able to use the 


4, | Next step is creating our Solid Geometry model. 
Click again on New Entity > Generation > 
Extrusion, 


In the construction by extrusion window, fill in the 
gaps according to the figure on the right. 


Your geometry also should have a look like in the 
figure, 


Finish with “Apply and Close” 


7 Res arm 
neme[Feed 
[Face creation from wites and/or edges 


objets @[[crees 


ty to create a planar face 


[Ears] arpy 


lose Help | 


lew ce 
SE 
the [Eset 
‘Base Shapes + vectxe ——— 
exe |= 
wie Le] 
iene: [F099 


au orectons 


F Rovere 
FF Seale the face oppoot tothe bss 
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EDDITING_COMMAND_FILE 


# |Description [Figure ‘I 
5. | We now have created our geometry and the following step | | ~ Create Group x 
would be defining our groups that will help us make the shape Te 
FE Analysis. oe ce ea 
Make sure Extrusion_1 is selected in the Object Browser. 
Click on New Entity > Group > Create en rere 
In the Create Group window that pops-up select the third ae. 
option for shape type and select the faces from the i 
geometry. : === 
This first group will be used for support, so it will help if oo oo 
you name your groups, and general all your options in the || | | M8” St#PeSeeten restriction Show only selected 
Object Browser, with easy to remember names. idee Hide selected 
L used “Fix” as my fixed support group. Oi | 
Fill in the gaps according to the figure on the right. | © Only Sub-hapes ofthe Second Shape | —O* SESE] 
Rotate the model so you can find the face where we need secons shape <@|[ 
our fixed support. When you click on the face we are 
interested in, it will get highlighted with white edges, Select Al 
Clicking on the “Add” button, the face number will be MM 
added as shown on the figure, aie 
Repeat the same procedure for selecting the opposite face 
and name it “Load”, 
Finish with “Apply and Close” 
6. | Now your Object Browser should look like in the figure on the 
right. 
The eyes on each option hide/show those options in the viewer 
window. = 
2 
2 
a 7 OZ 
> 8 / Circle 
2 @ MM Facel 
>  &@ Extrusion_1 
* oz 
* Face 1 
z ea Fix 
é 
7._| Now that we are finished, save your work to a destination of your choise and name it accordingly, 
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Step 4: Meshing Geometry 


# Step by step example of how to mesh our geometry is shown below, taking in consideration that you have completed 
our geometry and have saved it to a logical path with a meaningful name 


# [Description Figure 
7. | Pull down the category box or click on the mesh icon to [7+ aaa 


nter to Mesh Module. Name Mesh_1 


Geomesry [@ [Ensen T 


Select Extrusion_I from the Object Browser and click 
Mesh > Create Mesh and fill the pop-up window 
likewise, 


3p | 2] w| oD} 


Algorithm [Retgen 102050] 
In the Hypothesis tab, you won't have to change Hypothesis [NETGEN3DParamares =] aa] 2 
‘thing, just keep the data default. 
‘Aa Hypothesis [Toe SSCS at 
‘Accign a sot of hypotheses 
Finish with “Apply and Close” || #opivendGose| Apply | __ cose | Helo 
2. | Right-click on Mesh_I > Compute. 7 een 
<a 
If everything goes well, you will have the same pop-up window which ||| (= 
includes the mesh information. Pa 
es90 1 a 
rages 8 
Click on “Close” ||| rivimiors 0 
Your mesh should look like the figure below 


Ey Cap rer_Il: CYLINDER_ROD EDDITING_COMMAND_FILE 


Step 5,6,7,8 and 9: Salome_Meca Linear Static Analysis Wizard 


# In the previous chapter, in these steps we chose to use Salome_Meca Linear Static Wizard to create the command 
(.comm) file which would hold information about the analysis we wanted to do. BUT Linear Static Wizard doesn’t 
give us the option to set our force in an horizontal way. In this case scenario we will use the Linear Static Wizard 
to create our command file and then we will modify it. 


# |Description Figure 

1. | Click on Aster Module and make sure Mesh_1 is selected in the Object Browser. 

In Salome_Meca main toolbar, click Aster > Wizards > Linear Static and then a window will pop-up. 
Select “3D” as a type of Analysis, continue with “Next”. 

Select Mesh_1 and choose “Use of geometrical groups”, continue with “Next”. 

Enter 2e05 (MPa) for Young’s Modulus and 0.3 for Poison’s Ratio, continue with “Next”. 

For Geometric Boundary Conditions select “Fix” and give 0,0,0 values, continue with “Next”. 

For Load Boundary Conditions, it doesn’t matter what value you input here because we are going to modify it in 
the next steps, continue with “Next”. 

Save the file by giving it an appropriate name. 

“Finish”. 


# Now we have completed the steps to generate a study and create the command (,comm) file as we did in the previous 
chapter, Our next step is to edit the command file and adjust our force correctly. 


Step 9: Manual editing of command (.comm) file 


# | Description Figure 
7. | Search and open the folder where you saved the command [J Ta Feta 
file, it should be in the same folder where your -hdf file is <a 


4+) wane es Dep cere tn | sabe es) 


located. 


If your Ubuntu installation is unchanged, this file loads 
with “medi” text editor, by default, 


me 
a Pens 


igen 
2, | When the text editor opens, it should look like the figure on the right. 


You don’t have to worry or get scared what all those lines after lines 
mean or how they interact with Salome_Meca or Code_Aster, yet. WATL-M0OL_AILLAGE( reuse “HAIL, 
‘nTe PEA 2D="F (CAO MA= La, .9 


The modification we need to make in this step are really simple. Pe eR ea 


oDeLTSATTON='3D"=).) 


(epg Cap rer_Il: CYLINDER_ROD 


EDDITING_COMMAND_FILE 


Description 


[Figure 


# 
3. 


Find the sentence with the words “PRES_REP” which is 
the pressure term that we entered on face “Load” with a 
value of “1” (default). 


Select all the lines that relate to that command and delete 
them by pressing Delete. 


DEBUTC) 


MA=DEFI_MATERTAU(ELAS=_F(E=200000.0, 
eae. 


MATL=WODT_NATLLAGE(reuse =NAIL, 
MALLLAGE=MAIL, 
ORIE_PEAU_30=_F(GROUP_MA="Load" 


-AFFE_MODELE(WAILLAGE=WAIL, 
AFFE=_F(TOUT="OUI" 

PHENONENE= ' MECANTQUE! 

MODELISATION="30" ,), 


-AFFE_MATERIAU(MALLLAGE=MAIL,, 
(TOUT "Our", 
MATER=MA, ).) 


MECA_STATTQUE (MODEL E=NODE , 
‘CHAN_MATER=NATE,, 
EXCIT=_F(CHARGE=CHAR, ), }; 


CALC_ELEW(reuse =RESU, 
MODEL E-NODE, 
CCHAM_MATER=WATE, 


Replace the above sentence with the one shown in the 
figure on the right. 


“FORCE_NODALE” means Nodal Force in French. 


Here we replaced PRES_REP which may also be a key 
word to apply loads BUT it didn’t give us the option to 


apply that force in the direction that we were interested in, 


500N 


‘The amount of Load that we want to be applied is 
and it should have a direction of Negative Y (-Y). 
FORCE_NODALE > Nodal Force. 

GROUP_NO > Ina specific group of Nodes (Load). 
FY > Direction. 


DEBUT); 


WA-DEFI_MATERIAU(ELAS=_F(E=200000.0, 
MU=0.3,),); 


atl 


LTRE_MATLLAGE (FORMA 


MED" 


NATL=WODT_MATLLAGE(reuse =WATL, 
MALLLAGE=MAZL, 
ORTE_PEAL_3D="F(GROUP_MA='Load*,),); 


=_F(TOUT="OUr* 
PHENOMENE="WECANIQUE™, 
MODEL TSATTON='30', ),); 


-FE_MATERTAU(MATLLAGE-MATL, 
‘AFFE=_F(TOUT=" OUI" 


-FE_CHAR_MECA( MODELE=HODE 


‘OPTION=("STGM_ELWO" , “STEQ_ELNO",), 


excr 


F(CHARGE=CHAR, ),); 


Save the file with our new modifications. 
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Ste; 


p10: Running the Analysis 


# 


Description Figure 


L 


‘Open Salome_Meca again and check that your -hdf file is open. Make sure that Aster Module is selected. 
Right-click to the “linear-static” branch in the Object Browser and click “Run”. 

‘The pop-up window where the white letters and numbers flow will appear. 

If everything goes well, you will be able to see “Post-Pro” Module appear in the Object Browser. 


Most of the times that you run an analysis the “Post-Pro” branch will not appear in the Object Browser. That 
happens when there is a code error/mistake. To fix those mistakes you will have to go back to the command file, 
make your changes and then switch back to Salome_Meca, right-click on “linear-static” again, choose “Edit” and 
instantly click “Ok”, without changing anything inside that window. This way the analysis will get updated with 
the changes you made in the command file. Right-click and “Run” to compute the new data, 

There is a way that points out where exactly in your code appears to be the error, but we will address to that in the 
next chapter. 


Save your work, 


Ste 


p 11: Post-Processing of the results 


# Since you are familiar already with the way that Post-Pro Module visualizes the results, we will continue together 
but without thorough explanation. 


Description Figure 


‘Activate "Posi-Pro” Module from the drop down list - Deford shape +x 
Deformed Shape | sealar ear |. input | 


Open “RESU__DEPL” branch and add “Deformed Shape”. ———— 


See Factor: [EE] 
agate coarng 


To make the deformation more prominent, change the scale factor to 
“0.3”. 


Finish with “Apply and Close” 


For this Deformed Shape, enable “Show Only” and change the sss Sooo 
Representation to “SurfaceFrame”. opty and clse||aseiy | cancel |__| 
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EDDITING_COMMAND_FILE 


Es 


Description 


Figure 


2. 
27i82ie-20 29390 

a 

4 

5. | Don’t mind the results that are shown in the last two figures, The purpose of this analysis was to show you how to 


manually edit the command (.comm) file. 
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Chapter_IIT: Child Swing 


In this chapter we will examine the 1D_Analysis of a child 
swing beam. 
Emphasizing in Nominal Stress. 


Step 1: Purpose of the FE Analysis/ Description of the problem 


‘The purpose of this chapter is to find out the deformation of the beam, which consist the child swing, and the 
displacement of its end point. Also, we will examine how the stress will look like, when a loading by a specific weight, 
ned by the force. 


Step 2: Input Values for the FE Analysis 


Beam Dimensions: 1= 1000mm; 


(ey Cuaprer_Il: BEAM 


Step 3: Model Geometry 


# Open Salome_Meca, s 


1D_ANALYSIS, NOMINAL_STRESS 


a new study, enable Geometry Module and save your study with an appropriate name, 


# | Description Figure | 
1. | Click New Entity > Basic > Point. x Fein Constraction x 
———EEE 
A pop-up window appears like the figure on the right. os cs co 7x cm 
J : Fesult name 
Select the I* Algorithm in Points. sone 
Enter “P/” as a name. exit e$ @———_} 
=| 
Enter (1000,0,0) for X.Y,Z coordinates. =| 
| 
ey | ome | | 
Finish with “Apply and Close" | 
2. | Now we need to create a line with 2 points (0,0) and PI. 
i 
Click New Entity > Basic > Line, Gea ce 
; a 
Select 1 Algorithm for line by two points. oo 
Click on the blue curved arrow near Point 1, to choose our first | 
point of the line and repeat the same for Point 2, as shown in the rei @( 
figure on the right. 7 al 
Finish with “Apply and Close aaa Cie iw 
3. | Next phase is to arrange our groups, like we did in the > Create Group * | 
previous chapter. f Shape pe 
e+ oe ce 
Click on New Entity > Group > Create. Sane 
‘ : r a 
In the pop-up window we will follow the same steps we did aol 
before, BUT this time we will select the 1* Algorithm for Se $$$ 
selecting the points, which is located at the beginning of our Man shape @ [fires 
line. [ Main Shape Selection restriction _ shaw only selected 
Continue with “Apply” ||| | @ woresvcton ——— 
© sometrical pats othe sacond shape | ————"'_} 
Do the same for the other end and name it “Load”. © only Sub-shapes ofthe Second shape | Stow alsub-shapes] 
second shape _@ | 
Continue with “Apply” 
BO tet 
Now to create a group for our main line, we will need to sa 
select the 2“ Algorithm for Line, and name it “Guide”. jaan 
Select the line in the viewer tab and click “Add”, 
‘Apply and close] only _| tose Hep 
Finish with “Apply and Close" 
4,_[ Save the file and give ita proper name. 


yy Cuaprer_II: BEAM 
Step 4: Meshing Geometry 


# After we are done with modeling our geometry it is tim 
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¢ to move to Mesh Module. Enable Mesh Module in 


Salome_Meca and start the procedure to generate the Mesh. 


# |Description 


Figure 


1. | Click Mesh > Generate Mesh. 
Select Line_1 for Geometry. 


Select “Wire Diseretisation” for Algorithm from the drop 
down list. 


Click on the gear button next to Hypothesis. 


> Create mesh +x 
Nome Mesh 
‘seometry [@ [ine x 


so| 20 w| oo} 


Abe — TA =| 
Se ee] 


Le | 


Fe 


6 ypatacts 


2. | In the list that will fall down, select "Wh. Segments”, Keep 
default name and change the number of segments to “10”. 


Click “Ok”. 


Finish with “Apply and Close” 


3. | Right-click on “Mesh_I” from the Object Browser and click 
“Compute”. 


If everything goes well you should have the same pop-up wit 
shown in the figure, 


Apply and Coxe] apply | __ close |Help 
us Hypothesis Construction + xT) 
{BB Number segments 
“arguments 
Nome [We Segments 
Numbers Segments [13] 
Type of distribution [Equdiatant dainbuton =] 
OK Cancel Help 
compute men 
°@ 
indow as | 
ra! trae uate 
0 ements: 
ages wo » a 
wanes a ° 
uenges: 0 ‘ ° 
aise is 
Saat ‘ o 
Heap rns: 9 
teaeaaee 
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# | Description Figure | 
4, | Right-click on “Mesh_I” again and choose “Create Groups From 
Geometry”. Mesh 2 | 
‘Then you have to select the groups we created in Geometry Module com [ef 
and place them in the same box as in the figure on the right. cate 
— 
coonary 
Finish with “Apply and Close” | scene] soo = “sts 
# We are done meshing our geometry. Next step is creating the command file that will contain all the necessary 


information and will g 
our command file, 


e instructions to Aster about the analysis, This time we are going to use Eficas to conclude 


Step 5,6,7,8 and 9: Creating command (.comm) file with Eficas 


# As we noticed from the previous chapter there is more than one way to create a valid command file. In the previous 
chapter we used Salome_Meca Wizard and then we modified it in a text editor. In this chapter we will use Eficas 
Application. 

# Eficas is an application that generates syntactically correct command files, this does not mean that they will make 
sense or that there will not be any errors while running the analysis. This applications only makes certain that you 
won't miss any commas or any closing brackets, all other features have to come from the user. 

# [Description Figure 

T. | In Salome_Meca select “Eficas” trom the drop down list. 
A huge discomfort will be that eficas has not yet been translated to English and everything, including the tool tips 
are in French language. 

2, [In the main Toolbar, click on Eficas > Eficas > Tie arives x 
Icons from left to right: jo ee Siax 


1] Create new file, 2] Open file, 3] Save file, 4] 
Copy, 5] Paste, 6] Delete, 7] Cut. 


] 
Multicatalogue. i Elton Cte 
Click on Create a new file. | 
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Description [Figure 


# 
3. 


In the first pop-up window, you will be asked to : choieu code x 
choose a code. 


euler chois un code : 


Select “ASTER”. © Openturns Study 
© Openturns trapper 


siti 


Continue with “Ok” Cancel 


Tn the second pop-up window, you willbe asked to_| [ ~ Choir dune version du code ASTER x 
choose in which version of Code_Aster you want to 


3 versions du catalogue sont disponibles 
work on. 


Select “STAI” which is the newest one (at my 
system) and the most stable. Veuillez choisir celle avec laquelle zl 
Vous souhaitez travailler i 


ove | 


Continue with “Ok” | 


Before going any further, let’s save the file in the same folder we saved the study before. Click on Fishier > En 
registrer sous (Save us) and click “Save”. 


Tn the figure below you can see how Eficas application looks before we start working in it. 
In the left column it shows the file name “Coml.comm” in red colour which means that the file is not yet 
complete. 
Eficas is waiting for us to input the data, when that happens correctly the red square will turn green, meaning that 
our command file is correct. 
In the second column there are all the commands that are available in Code_Aster and are sorted Alphabetically. 
‘The data that needs to be added now, is called concepts, and for a fully functionable command file we need to add 
the following concepts: 
i, DEBUT [Start] 
ii, Material Definition [...] 
iii, What type of Mesh to read [. 
iv. Type of modelling [1D,2D,3D, Mechanical] 
v. _ Assign previously defined Material to the model [...] 
vi. Add Geometric Boundary Conditions [...] 
vii. Add Load Boundary Conditions [...] 
viii. Define Type of Analysis [...] 
ix. Results to be calculated at the elements, nodes [...] 
x, Save the results toa MED file [...] 
xi,__Fin [End] 
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# |Description Figure 
7. | Ajouter Mot-Clef = Add Keyword Mots Clefs Permis = Permitted Keywords 
Nommer Concept = Name of Concept Regles = Rules 
Nouvelle Commande = New Command 
Note: If you are having trouble with undestanding the meaning of some words/commands remember to use 
translation from French. 
& 
Ajouter Mot-Clet | Nommer Concept | Nouvelle Commande 
Mots Clefs Permis Regles 
9. | Find SDEBUTcommand from the second column, double-click on it to begin our command file. 
‘You can see that “DEBUT” square is green, which mean that this command is functionable and we can move to 
the next command. 
Clicking on the command that you just added, will get you back to the first screen where you can select your next 
type of command/concept. 
70. | Second step will be defining the material that we 
are interested in using. 
Choose “Nouvelle Commande” tab, to select anew | ia 
concept, and find “DEFI_MATERIAU” which a 
means DEFINE_MATERIAL. Double-click or (mee | 
select it and click “Valider” ee 
Now you can see that DEFI_MATERIAU has been as 
added to our command list but the square is yellow, aan 
which means that we need to add more information 
to make this command functionable. 
TI. | Find “ELAS” from “Mots Clefs Permis” and 
double-click on it to add it to the command 
file. enim 
awe other bjaneted - appear, 7 them like |; pea 
in the figure on the right, “E” stands for _— bree) 
Young’s Modulus and “NU” for Poison’s (eons = suc ceen 
Ratio. Comet 
After this step you will notice in the command on doneat: 
file list that every square is green except a 
DEFI_MATERIAU, and that happens because 
we need to assign a name to the concept from 
“Nommer Concept” tab. | 
72, | In the next step we must define the type of the | [ comtcomm | 
mesh that will be read as input. 
Commande Concept/Valeur 
Go to “Nouvelle Command” again and choose || ¥ @ coml.comm 
“LIRE_MAILLAGE” > FORMAT=MED and @ DEBUT: 
lastly give the command a proper name. Vv @ DEFL_MATERIAU : material 
> @ ELAS: 
y @ LURE MAILLAGE: mesh 
© FORMAT: MED 
® b format_med 
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#_| Description Figure ‘| 
73. | We came to the point where we will need to] [> 
apply Boundary Conditions on the nodes, but 2 
we can’t do that unless the MED file has Commande Conceptvaleur 
information for the node group. ¥ & comicomm 
@ DEBUT: 
The DEFI_GROUP command will help us do ¥ @ DEFL_MATERIAU : material 
that. > @ ELAS: 
v @ LIRE_MAILLAGE: mesh 
Click on Nouvelle Commande > © FORMAT: MED 
DEFI_GROUP > MAILLAGE > mesh. © bformat_med 
Click on DEFI_GROUP and double-click vy Gi DEFLGROUP: 
CREA_GROUP_MA, click on NOM and @ MAILLAGE: mesh 
name it “TOUT” and click “Valider”. ¥ @ CREA GROUP_MA: 
Click on TOUT > OUI which means give the @ Nom: Tout 
name (NOM) TOUT (all) to all mesh elements. © Tour: our 
¥ @ CREA GROUP_NO 
Here we made a group of the entire mesh and ¥ @ b option 
called it “TOUT”. OUT GROUP 0 
14. | Now we need to convert this group of mesh to a group of comt.comm | 
nodes. Commande Concept/Valeur 
BE com corm 
Click DEFI_GROUP > CREA_GROUP_NO > b_option DEBUT: 
> TOUT_GROUP_MA > OUI, which will convert all DEE MATERIAU maradal 
mesh group to Node groups. Behe a 
© FORMAT: MED 
@ b format_med 
y @ DEF GROUP mesh 
‘© MAILLAGE mesh 
¥ © CREA. GROUP_MA: 
© Nom Tout 
© tour: out 
¥ @ CREA GROUP_NO 
¥ @ b option 
© TOUT.GROUF_MA out 
Name the concept and move to the next one. 
15. | For the next step we need to define the type of modelling. ¥ i DEF GROUP: imesh 
“AFFE_MODELE” command is the most suitable Piaeeeee beam 
> @ CREA GROUP_MA: 
command for that task. §.@ CREA GROUP-NO. 
Choose AFFE_MODELE > MAILLAGE > mesh. 
Select TOUT > OUI, to asign the modeling to the whole mes 
mesh. oe 
Select PHENOMENE > MECANIQUE MECANIQUE 
Select MODELISATION > POU_D_T ¥ @ b.mecsnique 
Close the command with naming it. Sue eee Use 
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Description 


Next step will be to Add Material Properties 
which were defined previously to the model. 
Choose AFFE_MATERIAU (Assign Material) 
> MAILLAGE > mesh. 
Click AFFE > TOUT > OUI 
MATER > material (depends of the given 
name). 

Finish with naming the concept. 


17. 


Next step is to define the Geometry of the 
beam. 


Choose AFFE_CARA_ELEM and fill the gaps 
according to the figure on the right. 


Here we defined a Circle Cross-Section to our 
beam, with R=15mm. 


Finish with naming the concept. 


18. 


Tn the next step we will Add Geometric 
Boundary Conditions. 


Select AFFE_CHAR_MECA and fill in the 
necessary information as shown in the figure 
on the right. 


According to our data, we want 0 value to all 
degrees of freedom. 


Fi 


with naming the concept. 


19, 


In the next step we will Add Load Boundary 

Conditions. 

Select AFFE_CHAR_MECA again and give in 

the necessary information as shown in the 

figure on the right. 

Here we added a load of SOON to the “Load” 

group of nodes with negative FY direction. 
Finish with naming the concept. 


20. 


Next step is defining the type of analysis. 
Add MECA_STATIQUE (Static Mechanical 
Analysis) to the command file and then add 
information to the branch as shown in the 
figure on the right. 


Note that we are adding two EXCIT command 
to cater for two types of Boundary Conditions 
(EXCIT = excitation of the calculation). 


with naming the concept. 


[Figure a 
v GAFFE MODELE: model 
@ MAILLAGE: mesh 
> @ AFFE 
MAILLAGE : mesh 
our 
material 
v Gl AFFE_MATERIAU: mater | 
© MAILLAGE: mesh 
> @ AFFE: 
y G@ AFFE_CARA ELEM: Guide 
© MODELE model 
y @ POUTRE: 
© GROUP MA: (rour) 
© SECTION: CERCLE 
Vv ® brercle 
v_ @ bconstant 
CARA R 
© VALE: (15) 
v Gl AFFE_CARA ELEM: Guide | 
© mopeL model 
> @ POUTRE: 
| ARFECHARMECAS Five 
‘© MODELE: model 
v @ DDLIMPO: 
© GROUP_NO: (Fin) 
© px: ° 
0 
° 
° 
) 
° 
Fixed 
model 
Load 
MODELE : model 
v @ FORCE NODALE: 
© GROUP_NO: Load 
on: -500 
y GAFFE CHAR MECA: Load ral 
@ MODELE: model | 
> @ FORCE_NODALE : 
| MECASTATIQUES result | 
MODELE: model 
© CHAM_MATER: mater 
© CARA ELEM: Guide 
v @ exci: 
y @ exci: 
© CHARGE: Fixed 
vo exc: 
© CHARGE: Load 
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#_| Description Figure 
21, | Next we define the types of result to be etm eGR: =—* 
calculated from the initial displacement : 
aiid: @ MoDELE: model 
‘Add CALC_CHAMP to the command file and : Goren par 
fill the information according to the figure on pata 
None Je CALCCHAMP TT resutt | 
Name it with the same name you used with the Cer a 
previous concept. - pte 
Here we want results for 2 information i 
Calculated: af 5 arpa? ‘SIPO_NOEU) 
SIPO_NOEU = Stress at Nodes. ‘Sues: ; Caan 
REAC_NODA = Reaction at Nodes. z | 
22. | In this step we need to define the information 
that will be written in the MED file. MED 
@ b-format_med 
‘Add IMPR_RESU and follow the way the © b-restreint 
figure shows. Vv @ RESU: 
© MAILLAGE: mesh 
At this point our goal is to print and visualize @ RESULTAT: result 
some results. @ binfo_med 
® b partie 
DEPL for deflection and the other 2 that we © buvari.et 
explained in the previous step #21. @ bextrac 
‘@ bmp 
@ bmed 
@ btopologie 
23. | Last step will be adding the FIN command. coml.comm 
Once you add FIN, “coml.comm” square turns | | Commande _Concept/valeur 
green, meaning that our command file is ¥ Bi comi.comm 
completed. @ peeuT: 
> @ DEFI_MATERIAU: material 
> @ LIRE MAILLAGE: mesh 
> @ DEFILGROUP: mesh 
> @ AFFE_MODELE: model 
> @ AFFE_MATERIAU : mater 
> @ AFFE_CARA ELEM: Guide 
> I AFFE_CHAR MECA: Fixed 
> I AFFE_CHAR MECA: Load 
> I MECA STATIQUE : result 
> Bl CALC_CHAMP: result 
> @ IMPR_RESU: 


27._| Save this file and return to Salome_Meca. 
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#_| Description 


7._| Now that we created our command file with the 
help of Eficas Application, we need to define a 
study case in Aster that will do the analysis and 
point the location of our command file. 


Choose Aster Module from the drop down list. 
Click Aster > Add Study Case. 


Chose a correct name. 


Locate your command (.comm) file from your 
system, 


Set Mesh_I to the Mesh. 
And input the Solver Prameters as shown. 


Finish with “Ok”. 


| Figure 
cessation el 
Sy exe defi 
Nae [Smt 


ammandsle [Far dex =] Ja fore arms eet comer 


Mean [From object orowser zl @ | [mena 
Ast seviee 
ve abet _5f eee = of ee 


Execution made [iraaet? =] Interactive folaw up FF 


Solver prameters 


2._| When you click “Ok”, a new branch appears to the Object Browser 
named Aster, which contains the study case that we just created. 


Choose and right-click to “Run” the analysis. 


the figure on the right. 


= 
If everything is correct, after the analysis is complete you should be Be 
able to see “Post-Pro” branch added to the Object Browser, like in tine. 


Seen te) [Zee 
Bee sean 

coe | 

are ol 


& & Applied hypotheses 
1 & Applied algonithrns 


aa 
8 


b+ Beamismed 
& mesh 

Families 

@ Groups 

Rags 
 result_DEPL, 
E result_REAC NODA, 
 result_SIPO NOEU, 


3._| Save the study and let's have a look at the results. 
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Step 11: Post Processing the results 


see what the 


s has to show. 


# Open Post-Pro Module from the drop down list and let Analys 


# |Description Figure 

1. _| Follow the same steps we took in the previous chapter. 

Result_DEPL > Deformed Shape > Show Only > Surface. 

We can see that the deformation has a range from 2.32934e-21 mm to 20,9694 mm. 


resto 


3. | Next we will need to see the stress results, 
Result_SIPO_NOEU > Scalar Map > Mode: SMFZ (Stress due to bending Moment MZ) > Show Only. 
We can see that the direct stress has a range from -7.02694e-13 MPa to 188.628 MPa. 


resut_s0NOEUG, 


5._| Below you can see the figure which includes our command file with all of our concepts/parameters/data. 


6 | Dee resul t=-WECA_STATIQUE(OOELE=nodel | 
wnaterial-DEFL_MATERIAU(ELAS™_F(E-2005, cnRa_ELE 
NU=0.3..95 EXCIT=(_F(CHARGE"Fixed, ) 


FCCHARGELoad,},),) 


Ipesh=LIRE_MATLLAGE(FORUAT= "WED", ) 
result=caLe chan rus 

‘CRITERES="S1EQ_NGEU', 

FORCE="REAC_ NOUR") 


ipesh-DEFT_GROUP(reuse = 
MATLLAGE=a0sh, 
(CREA_GROUP_A_F(NOW= "TOUT" 
TouT="Our";) 
(CREA_GROUP_NO=_F(TOUT_GROUP MA"OUT", ),) 


1MPR_RESU(FORUAT= "NED, 
RESUS F(MALLLAGE*nash 
RESULTATeresult, 
NOW_CHAU=(“DEPL', "REAC NODA" 'SIEQ_NOEU',“SIPO_NDEU'.).9,) 


‘node “AFFE_NODELE(HATLLAGE-nesh 
‘AFFE=_FCTOUT= "OUT 
PRENEMENE="EcANQUE 
MoDELTSATION="POU_D-7",),); 


FIN) 


imater“AFFE_NATERIAU(UATLLAGE=nesh #CHECKSUM:1250ebeBacca44gd92118826032e89¢2 -:FIN CHECKSUM 


ATER=nater al, ), 


[Fixed=AFFE_CHAR_MECA(MOOEL 3 
DOL_IMPO=F(GROUP_NO= Fix’, 
02-0, 
R=0. 
Dav-0, 
DRz=0, > 


\FFE_CHAR ECA(WODELE-node).| 
FORCE_MODALE™_F(GROUP_NO-'Load 
FY=-500,),) 


Ey CaP rer_IL: BEAM 1D_ANALYSIS, NOMINAL_STRESS 
Notes 


i, Code_Aster can read IDEAS, GMSH, MED Formats. 
If in Eficas you are having trouble understanding a word, use translation from French to English, it is really 
helpful. 

iii, In defining material step, we filled only 2 parameters (Young’s Modulus and Poison’s Ration) because the task 
didn’t require extra ones like Density, Thermal Conductivity, ete. 

iv. There will be cases where you will have two different material types. To do so you will just create another 
DEFI_MATERIAU concept with different parameters and you will just give it another name. 

v. This can happen to other commands also, like in AFFE_MODELE. But have in mind that you can separate groups 
and give different parameters to each one of them only by changing their names. 

vi. The reason why we define mesh group is that we can convert it to node group later. 

vii. We need to be careful of which type of modelization we choose, because in most of the cases we pick the wrong 
one and the command file is not working. To choose the best one in a certain case study, you will have to study 
carefully the Code_Aster documentation, U4.41.01. 

viii. _ Beams and Shells have 6 Degrees of Freedom, 3 in Translation and another 3 in Rotation. 

ix. In Aster, regarding the default 128MB memory parameter, you could leave it like this but soon, for bigger 
analysis, that amount of memory will be insufficient and it will stop running. The more memory you input the 
safer it is. Have in mind that the most sufficient amount of memory is around quarter or half of the total memory 
of your system, 

x. There will be times that your analysis will not allow you to visualize Post-Pro Module, that is happening because 
there is an error in your command file. 

xi, When you Run an analysis, it does not matter if it is working or not, Salome_Meca creates some files in the 
folders directory. [.hdf, .rmed, .resu, mess, .mmed, .export, .comm file]. 

xii, In the .mess file as we mentioned before, there will be information regarding your analysis. At the end of the file 
there will be a “box”, in which will be stated the area, the line where your code is mistaken, In the next chapters 
you will see examples of that file. 
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Chapter IV: Beam System 


In this Chapter we will deepen more into 
beams, 

We will study a 1D be 

emphasizing in Nominal Stress and Alignment. 


m system analysis, 


Square Cross-Section 


Step 1: Purpose of the FE Analysis/ Description of the problem 


In this chapter we are asked to create a 3D Cube Beam System. The Cube Beam System has a square Cross- 
Section which is given above 


Step 2: Input Values for the FE Anal 


Cube Dimension: 100mm x 100mm x 100mm; 
Cross-Section: Rectangular, Length= 8mm x 8mm; 
Young’s Modulus: 2e05 MPa; 

Poisson’s Ratio: 0.3: 

Pressure Applied: Fl=150 


200 


curs 
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Step 3: Model Geometry 


# 
# 
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As we proceed forward and you learn more and more things, the explanation of the steps will be shorter. 


Start Salome_Meca and begin modelling. 


Description 


Starting with, we must add all cube parts in Geometry Module. 
Start Geometry Module, click on New Entity > Basic 
Point. 


Enter the following Points: 
Name = Coordinates (x, y, 2) 

P1_0,0,0 = (0,0,0) P5_0,0,100 = (0,0,100) 

P2_0,100,0 = (0,100,0) P6_0,100,100 = (0,100,100) 

P3_100,100,0 = (100,100,0) P7_100,100,100 = (100,100,100) 

P4_100,0,0 = (100,0,0) P8_100,0,100 = (100,0,100) 


Finish with “Apply and Close” 


Figure 
. Point Construction x 
(Points 

fe Ce Fre ce Coe 
- Result name 


Name [FI_0.0.0 

/ Coordinates 

CS | 
yf | 
2 [or | 


[marie] ay | cme |_ mew | 


‘After creating the Points we must make Lines out of those 
Points. 


Go to New Entity > Basic > Lines 


Connect all Points according to the figure on the right and our 
given data. 

Line_13 = (P5,P7) 
Line_14 = (P8,P3) 
Line_15 = (P4,P2) 
Line_16 = (P1,P6) 
Line_17 = (P6,P3) 


Line 
oe Ce 
a 
Nome [lined 
a 
voit @| [Bae 
roit2 [je [203000 
apply | Chose | Heb 


After finishing with the Points, Lines you should be able to see 
the geometry as shown in the figure on the right. 


‘Object Browser ax 


TE @® Geometry 
ro 


P3/100,100.0 


CO) 
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Description Figure 
Next we need to create a container that will store all these 
individual lines. -compauna 

°@ 
Click New Entity > Build > Compound. as 

fame [componnd. 
Keep default name. cone 

Sasa 


Select all our Lines and Points for Objects. 


Finish with “Apply and Close” 


— a 
nt oat 


TT Set presentation parameters and sub-shapes from arguments 
F Add preicto names of restored sub-shapes 


Next we need to create Groups of Nodes for our Geometric and 
Load Boundary Conditions. 


Click New Entity > Groups > Create. 


In the window that pops-up, select the 1" algorithm for Shape 
Type of Nodes. 


Enter “Fix” as a name (or any name that suits you) and select 
the Points where we want our fixed support to be placed, from 
the Object Browser. Then click on “Add” and the Points will be 
added likewise the figure on the right, 


Click “Apply” 


Similarly create Node group of our loads, for Load_F1 and 
Load_F2. 


Click “Apply” 


7 reste Group x 
Shape ype 

a c co ca 

‘Group ame 


Name [Fa 


Main Shape And Sub shapes <<< 


Mein shape _@ [[Cempound 


Select now the 2 Algorithm for Shape Type to select Edge 


Name it “Guide” (or any name that suits 
All” and then on “Add”. 


you), click “Select 


All Lines/Edges will be added now on the list like the figure on 
the right. 


Finish with “Apply and Clos 


-Main shape Selection restriction Show only selected 
No rerction paar eaeaa 
ide celeste 
(© Geometreal parts ofthe Second shape 
© Only Sub-shapes of the Second Shape SNOW al sub-shapes] 
Second Shape @ =| 
Select al 
id 
‘apy and Close] Asply He 
: croup x 
shave pe 
a G ca re 
‘rup ame 
name [oude 
Man Shape And Sub-shapes 
Main shape _@ | [Eonpon 


Hain Shape Selection resection 


=a 
pares 

as 
eS eee |e 
© only Subshapes ofthe second shape Stow allsubyshapes] 
seo —$j 

aa 

7 

(everce] ery | ome | we 


This concludes our Geometry creation session. Save in a proper location with a proper name. 
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Step 4: Meshing Geometry 
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#_| Description Figure 
T._| Select and open Mesh Module from the drop down list. : creates ts 
Click Mesh > Create Mesh (with Compound_t selected). Name Weak 
a ~~ 
Select “Wire Diseretisation” for 1D Algorithm. 
| aw} ol 
Select Nb. Segments_1 with number of Segments “/0” as Agcrshm — [itreDecrtenion Sd 
Hypothesis. twypathese [ih Seqmens1 5] af 
a, Hyptvesis [Emme 
sign a eo f 
Finish with “Apply and Close” | | Aesivsedess|_—_toviy _| |___ ee 
2._| Right-click on Mesh_1 > Compute. Ks at 
ore 
If everything is correct you will have the same results as in the | | © 
figure on the right. Res 
Mesh 
We know that we have 17 edges/lines with 10 Nb.Segments, ||" on 
so there should be 170 total edges. in 
Fanaie ° ° ° 
aoe b ‘ ° 
ieee ane 
Se a 0 
Eee 6 ° ° 
feabeioes:  o a ° 
See ° ° 
pene ° ‘ i 
tesapoa pri © 
eee ie 
Click “Close” = 
3._| Right-click on Mesh_1 > Create Groups From Geometry. 
Select “Fix”, “Load_F1”, “Load_F2”, “Guide” from 
Geometry’s Object Browser and the groups will be added 
likewise the figure on the right. 
F Nodes 
caomy €| 
Finish with “Apply and Close” | topww=ra coe] __sp0v_|__Gose te __| 
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Step 5,6,7,8 and 9: Creating the command (.comm) file 


# Up until now we have seen a couple of ways to create the command file. At the beginning we used Linear Elastic 
Wizard from Salome_Meca, later we used the help of Eficas Application, Now we will edit by hand in a text editor 
an existing command file (the one from the previous chapter), and we will modify it to fit our new demands. 


First thing you need to do is Copy the command file from the previous chapter and Paste it to the folder of this case. 
Rename it with an appropriate name and open it with a text editor. 
Red letters are the new modifications and 


ae 3b ak 


Are commentary to help you understand what we are doing. 


1. DEBUT (); 
2 


3, material = DEFI_MATERIAU ( ELAS= _F ( E = 2e05, 


4, Nl 

5, mesh = LIRE_MAILLAGE ( FORMAT ='MED, ); 

6. 

7. mesh = DEFI_GROUP ( reuse = mesh, 

8. MAILLAGE = mesh, 

9. (CREA_GROUP_MA =_F (NOM ='TOUT, 

10. TOUT ='OUI, ), 

I. CREA_GROUP_NO=_F ( TOUT_GROUP_MA ='OUT, ), ): 


12, 
13, model = AFFE_MODELE ( MAILLAGE = mesh, 


14, AFFE =_F (TOUT ='OUI, 

1s, PHENOMENE = 'MECHANIQUE,, 
16. MODELISATION = 'POU_D_T’, ), ); 
17. 

18, mater = AFFE_MATERIAU ( MAILLAGE = mesh, 

19, AFFE = _F ( TOUT ='OUI. 

20. MATER = material, ), ); 


21, #We now want rectangular Cross-Section, so we need to change CERCLE --> RECTANGULAR and define our new 
22, dimensions. 


23. 

24, Guide = AFFE_CARA_ELEM ( MODEL = model, 

25. POUTRE = _F (GROUP_MA = 'TOUT’, 
26. SECTION = 'RECTANGLE’, 
27, CARA = ( ‘HY’, 'HZ 

28. VALE # VALE = Values 
39. 

30, Fixed = AFFE_CHAR_MECA ( MODEL = model, 

31. DDL_IMPO = _F ( GROUP_NO = Fix’, 
32. 

33. 

34. 

35. 

36, 

37. 


38. 
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39, #In this case study we have different load cases. We got 2 different forces now Fl and F2 (Load_F1/F2). 
40, #Keep FORCE_NODALE concept but edit the parameters. Load_FI in -X direction and Load_F2 in Y direction, 
41. #Magnitude F1=150 N, F2=200 N. 


42. 
43, Load = AFFE_CHAR_MECA ( MODEL = model, 

44, FORCE_NODALE =(_F (GROUP_NO = 'Load_FI’, 
45. FX =-150,), 

46. _F (GROUP_NO = 'Load_F2’, 

47. FY = 200, ),),) 

48. 

49, result = MECA_STATIQUE ( MODEL = model, 

50. CHAM_MATER = mater, 

Sl. CARA_ELEM = Guide, 

52, EXCIT = (_F (CHARGE = Fixed, ), 

53. _F (CHARGE = Load, ), ), ); 

54, 

55, result = CALC_CHAMP (reuse = result, 

56. RESULTAT = result, 

57. CONTRAINTE = ’SIPO_NOEU’, 

58. FORCE = ’REAC_NODA’, ); 

59. 

60. IMPR_RESU = (FORMAT = 'MED’, 

61. RESU = _F ( MAILLAGE = mesh, 

62. RESULTAT = result, 

63. NOM_CHAM =( ‘DEPL’, ‘REAC_NODA’, ’SIPO_NOEU’, ), ), ); 
64. 

65. FIN (); 


# Save this file and run the FE Analysis now that we have input all the required data. 
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CHAPTER_T 
Step 10: Run the Analysis 


and a command file. Since we have both of those 


# Atthis point in order to run the FE Analysis we need a Mesh file, 
files let’s start the analysis and then proceed to the visualization mode. 


# |Description Figure 
1. | Activate Aster Module, go to Aster > Add Study Case |= bape 
Fill in the gaps according to the figure on the right shay co ein: 
Give it a considerate name, choose the command file 
from its folder location, choose the mesh file from the epee oonos 2! | 
Object Browser 
ten romance of |_| [pen 
Change to STALI of Aster Version, and lastly increase 


the Total Memory (MB) and Time (s). 
a Teeatom =] Anerverion  [STRIT —=] tea svar 


execution mode [Free] interactive flow tp 
‘She preter 

Zeralmerony ad [aces = ine a a0 
cru number [TS ve esa trance 


cance 


Finish with “Ok” 
2. | Anew branch named “BeamSystem” (the choosen name) will be created in the Object Browser 

Right-click on it > Run, to start the Analysis. When the analysis is done and everything has been done correctly 
a new branch named “Post-Pro” may appear. 


Step 11: Post Processing the results 


# Now that the FE Analysis is complete let’s have a look at the visualized results Post-Pro Module is offering us 


# |Description Figure 
We will not do anything new in this step, that we haven't already done. 

ctivate Post-Pro Module from the drop down list, open result__DEPL branch, right-click on 0,- > Deformed 
Shape > tick Magnitude coloring > “Apply and Close” > right-click on Defor.Shape > Show Only > 
Representation > Surface. 
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# |Description Figure 
3. | Same process for the Stresses also. 

Open result__SIPO_NOEU > Right-click on 0,- > Scalar Map. 

In this case we will u Scalar Mode. 

Right-click > Show Only. 
4 


103742 


ey Cuaprer_V: PLate 2D_ANALYSIS, PLANE_STRESS, NOTCH, SINGULARITY 


Chapter_V: Plate 


In this Chapter w 


vith a notch in the 


mcr 


Step 1: Purpose of the FE Analysis/ Description of the problem 


_V: PLATE 2D_ANALYSIS, PLANE_STRESS, NOTCH, SINGULARITY 


In this chapter we are interested to know the concepts of Finite Element Simulation, like “Stress Discontinuity” 
“Structural Error”, “Stress Singularity”, “Finite Element Convergence”, To demonstrate this concepts we are going to use 
a Plate geometry with notch, which is made out of steel and has the given dimensions. The bar has Fixed Support on the 
left side and is subject to a tension of SOMPa on the right side. 


Step 2: Input Values for the FE Analysis 


Dimension: Length = 100 mm [1]; Length = 25 mm [b]; 
Length = 50 mm [a}; Thickness = 10 mm [t]; 
Radius = 15 mm [R]; 


Young’s Modulus: 200 GPa; 
Poisson’s Rati 03; 
Force Applied: 50.N[F}; 


comcar 
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Step 3: Model Geometry 


PLANE_STRESS, NOTCH, SINGULARITY 


# Open up Salome_Meca and start the input for Bar Shell Model, 
# |Description Figure | 
Z._| To generate a Shell model you just need to create a planar face. = 2 Sah Coats = 

connate sytem 
Click New Entity > Basic > 2D_Sketch, [Saba comtnare apron =] _ restore 
¢| =n memos 
In the window that pops-up keep “Global coordinate system”, 
“Element Type” and “Point” default. 
ony ca 
Enter for first value, Pi) > X=0_ Y=0 
Click “Apply” 
Your next points should be: P2 > X: reine 
P> 6 Able 
Bi s 6 rele 
Ps > X=200 _ Y=25 € section 
P, > X=200_ Y=0 
aes 
“2 a 
a cd | ere 
Click “Sketch Closure” 
gece | starch cowure] cancel | _ He 
2. | We need to create a mirror image. ~ rein Obra «| 
ero 
Go to Operations > Transformation > Mirror Image. ce ee ct 
7 ——_————SSS= 
Select 2 Algorithm for Mirror. a 
gues 
ones @|[eana 
Axe Mirror: OX Plane. coteces [EI 
FF crete scopy 
Select “Create a copy”. vanced opine 
1 Set presentation parameters ac sb shapes ra aguments 
FF a pret nares of strd sub-shapes 
review 
Click “Apply and Close" | |_ === aaa ea 
3._ | Next you need to convert this Wires into a Face. > ‘reate A Face ail 
rice 
Go to New Entity > Build > F: ca 
Result nam 
Keep default name, and in the Object Browser, select Sketch_1. a 
emefeet 
Click “Apply” | | =Fseeereston from wines andloredges———— 
Do the same for Face_2 with Mirror_1 for “Objects”. objets (@ [San 
ry torte planar face 
Select “Try to create a planar face”. 
[feeysricece] sop | coe | neo 
Click “Apply and Close” (Eeenes eI Ee) 
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Description Figure ] 
Next part is to create a shape out of the two faces we assigned. | Fare Two ojecs x 
ruse 
Operations > Boolean > Fuse will do the trick. a 
. — 
Leave default name, and assign the Arguments accordingly. a 
Face_1 for “Object 1”. Arguments 
obet1 _@ | [RE 
Face_2 for “Object 2”. | 
‘advanced options 
TF Set presentation parameters and sub-shapes fram arguments 
FF dd prefs to names of restored sub-shapes 
Finish with “Apply and Close” | | [Besivané es] _ son'v_| cose | rer | 
Next part will be, creating our Geometry groups, which | | > Create Group x 
will help us to assign supports, load cases, forces, etc. shape Type: 
c- ce 
Go to New Entity > Group > Create. 
‘Group Name 
In the Create Group window that pops-up, select 2" fa 
Algorithm for lines, eo 
" select in Shay [Fuset 
Enter Fix as Group Name, and then select and add each Min shape _@ | [Tuse2 
one of the two left lines where we want our Fix Support ons Shia Setar tetration Shaw only selected 
to be. We got two lines here and not only one because © Na resection ee 


we used Mirror Image to create our Geometry. 
Click “Apply”. 
Likewise select the opposite lines and name them Load. 


Change the Shape Type algorithm to the 3 option, the 
Face selection. 


Name your group Guide and Add the two face available 
to you. 


Finish with “Apply and Close” 


© Geometncal patts of the Second Shape 
© only Sub-shapes of the Second shape 


ome 


u Select All 


‘Show all sub-shapes 


Remove 


[ssi and cose] snow _| gee | _ tow _| 
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Step 4: Meshing the Geometry 


# In this step, we are going to use all of our resources to set up the mesh and then create our command file (comm). 
Let’s get started by opening Mesh module from the drop down li 


#_| Description Figure 
T._| Go to Mesh > Create Mesh. > ey 7 
Name [eh 
Fill in the gaps like the figure on the right. Geomeiy [@ [Ret 
But this time when it is time to input Max. Element Area_l, ao 20] 10 | oo] 
change the default value of Max. Area=20. — fa __——€| 
Hypxthese [fax Eanet arma =] af 
For ID tab, choose Wire Discretisation for Algorithm and oS 
Hypothesis would be Max. Size=10. ioMecis, BTES a ala) 
pilin sign a sat ofhypoteses i 
Click “Ok” | | ee 
Finish with “Apply and Close”” 
2. | Now right-click on Mesh_1 > Compute. Mesh information. PL 
: vo ier ane 
window will pop-up. eis aa 
Right-click on Mesh_1 and select “Convert from/to quadratic”. | | euem™m=i_—® 
Click “Apply and Close” || 8 - = 
angles: ° na 
Compute Mesh_1 again, and you will have results according to | | ousionale: 0 ° ° 
the figure on the right. robaene: 8 i 
3._| Last thing we need to do, is to create groups from our Geometry. [7 > create croup from ceemety +x 
- Mesh Le Pie 
Right-click on Mesh_l > Create Groups From Geometry. ae 
Select “Fix”, “Load”, “Guide” from the Object Browser and add Ea 
them to the window as shown in the figure on the right. 
Unselect the blue curved arrow next to “Geometry”, to be able to 
select Fix for the Nodes. 
ihe 
Geometry la 
Finish with “Apply and Close” | #2304 92s] 200) au = 


4_| Save your work, and continue to the next step. 
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Step 5, 6, 7, 8, 9, 10: Creating the command file, run the analysis 


# 


Although we gained some experience in creating the command file with a couple of different tools, opening a text 


editor and start coding, is s 


ill a long way down the road, Eficas Appl 


ion will do the trick again, and we will 


also use a text editor to edit some points, to get more familiar with coding. 


Description 


Figure 


‘Activate Eficas module from the drop down list. Go to Eficas > Eficas MultiCatalogue. In the pop-up window, 


click on Create New or Nouveau. 


Choose ASTER > “Ok”, STAI] > “Ok” and you have set up Eficas and it is ready for modification. 
Before starting editing, save the command file in the same location where the geometry/mesh/.hdf files are 


currently located, with an appropriate name. 
Fishier > Enregistrer sous. 


The figure on the right shows the first commands that 
we need to conclude in the command file. 


DEBUT > Begin. 


DEFI_MATERIAU > Define the material, Structural 
Steel is the material we will use in this case study. 


LIRE_MAILLAGE > Read the mesh, with MED 
Format. 


MODI_MAILAGE > Modify the mesh. 
ORIE_PEAU_2D > Reorientation of mesh, to be able 
to apply pressure on a “Skin”. 


PlateBarCamm.comm 


‘Commande 
¥ Gl PlateBarComm.comm 
@ peeuT 
y Gl DEFLMATERIAU 


Concept/Valeur 


material 


AFFE_MODELE > Apply/ Assign to the model, the mesh, 
phenomenon and type of modelization. In our case we are dealing 
with 2D plane elements that is why we use C_PLAN. Also we 


want Mechanical phenomenon. 
AFFE_MATERIAU > Apply/ Assign material. 


AFFE_CHAR_MECA > Apply/ Assign Charactei 


with Degrees of Freedom. In this case we want Fixed Support on 


the “Fix” group. 


AFFE_CHAR_MECA > Apply/ Assign Characteristics, Load 


cases. Pressure of 50 MPa to group “Load”. 


© MAILLAGE: mesh 
¥ @ ORIEPEAU20: 
© GROUP_MA: Load 
Fotesorcommeomm | 
Concepuvateur 
material 
mesh 
mash 
model 
mesh 
vO ASE 
© GROUP MA (Guide) 
© PHENOMENE: MECANIQUE 
¥ @ bimecanique 
@ MODELISATION: cPLAN 
WL AFFE.MATERIAU. materiau 
@ MAILLAGE: mesh 
¥ @ Are: 
© Tour out 
© aren: material 
17 AFFE_CHARLMECA: Fixed 
@ woDEE adel 
¥ @ DoLiMPO: 
© crour No Fa 
© ox. ° 
© oy. 0 
@ oz: o 
© ex o 
© pry ° 
© pez ° 
1 AFFE CHAR MECA: loade 
@ mopete model 
¥ @ PRES REP 
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# |Description Figure ‘| 
4, | MECA_STATIQUE > Linear static mechanic with FlsteBarComan. comm 
superposition on loads. commande oncepalour 
YB Patenarcomm comm 
CALC_CHAMP > Create a result by calculating fields by Sane iceaas ae 
elements or nodes. The “b_lineaire” branch is the command || + m unesatiace: mesh 
that allows Code_Aster to calculate specific results. In this Pa Monae mah 
. > Are MODELE mode 
case we are interested in Averaged (NOEU) Normal Stress || > m avte.warctiny materi 
(SIGM_) and Unaveraged (_ELNO). > AFFE CHAR MECA: Fae 
> mL AFE_CHAR MEA Loade 
McA stavique recat 
POST_RELEVE_T > This commands summarizes @ wooeie model 
assigned results and puts them to the .resu file, in a tabular = aria 
or a conventional format. In this case we are interested in 
summarizing the reaction forces. The first abstract of code, aed 
puts the sum of the individual reaction on the group “Fix”, ree 
ina table named “reakce”, with ym cuc.cHANe result 
OPERATION="EXTRACTION’ keyword. pAb eea ba 
{STON ELNO.STGMLNOEU) 
yl POST RELEVET reakce 
¥ @ ACTION 
© orenanow: extraction 
fe inmmue sum reactions 
ESULTAT: res 
¥ @ bexrac 
NOM. Cuan REAC NODA 
TOUT_ORORE ou 
$b pre.rla 
vo baure 
tmxoy 4 
5. | IMPR_TABLE > Prints the table of summarized reactions. | /m=aiGenan 
IMPR_RESU > Stores down and prints the final results. Potp= nee 
The “b_extrac” branch creates the branches in Post-Pro pee feo 
module. ootianir ie 
FIN > Finish. poneneners ery 
inca sie rest 
Bt rosea ice 
Stine ae 
bes 
hares 
ron eo 
$ Siomates 
ohne 
+S nas 
© b info_med ee 
eave 


(Der stP0_ENO ston EO To, NOEL) 


6._| Save the command file, 
In the next parts, we will check if the command file we created is functional, otherwise we are going to locate our 
mistakes and try to fix them. 
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# |Description Figure 
7. | Close Eficas Application and get back to . abppan = 
Salome_Meca. sy cso ttn 


Name fee 
Activate Aster Module from the drop down list, go to 
Aster > Add a Study Case and fill in the gaps oe Ee) [OT 
according to the figure on your right. 


Bat [rom opeccbowre =] rn ater 
Click “Ok” 
os 
A new branch appears in the Object Browser with the sever Taahos a] Atervesen [STAT I) ashame 
same name you indicated before. xen mode [rect] Intec fon uo 
SS 
Right-click on it and select “Run’ es mane (se) [TS] ne) — 


umimbe FT Sweet atanie 


om] cone | 


&. | As you noticed, the analysis is done but Post-Pro branch did not appear like it was supposed to. This is because 
there is a mistake or let’s say a wrong input in our command file that needs fixing. To do so we will have to find 
a file with the extension .mess, which contains all the analysis information including the errors, 
This file is located to the same folder with your .comm and .hdf file. 
After locating the .mess file, open it with a text editor and scroll your way down. At the bottom of the text, you 
can see a “box” created with exclamation marks, Inside that box lays the information about what went wrong 
during the analysis. 
One of the most disturbing issues you will face with Salome_Meca and Code_Aster is that everything is in 
French, including the analysis information text file (mess). 
9. | According to my mess file, there is an error in the J] 4. egecun a uexecution ~ rwreeRurTION 

DDL_IMP command, regarding DRY orientation, | |>> Joc.py + pebut farPorT 


es geet (Gk avesecution de JDC en MIXTE 
Itis fair error if we take into consideration that we | |“* 4/¢**eutien de 0c en WIXTE, 


are dealing with 2D Plane Elements, meaning that || 1 <s> exception utitisateur levee aais pas intercept 
DRY,DRX,DRZ and DZ orientations are Heap efron rss 
redudrant. ! ' 
1 aucun nocud ne connatt e opL: ORY 

fin CR d' execution de JOC en MEXTE 

>> J0C.py + FIN eaPpORT 

EXECUTION CODE_ASTER_EXIT_0006-5629-netwe-VirtualBOx=1 

70. | Open your command file with a text editor, and | | Fixed=AFFE_CHAR_MECA(MODELE=model, 

locate the part where your fixed support is. DDL_IMPO=_F(GROUP_NO=" Fix", 


Delete manually the hightlighted orientation as we 
will not need them in this Analysis. 


Loadz=AFFE_CHAR_MECA(MODELE=model , 
PRES_REP=_F(GROUP_NA='Load", 
PRES=50.0,),); 


71. | Save the changes in the text editor and get back to Salome_Meca. Now we need to run the new Analysis again. 
In order to update the stream data, right-click to your study case > Edit, leave everything as they are, right-click 
again > Run. 

What we did now, is update the study case with the new modification we made. 
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# 


Description [Figure 


2. 


‘Apparently we got more than one code anima 


mistakes. This time, if you translate the 

message on the right, it says that 6 mesh 1 ype de lvexception = 
points do not belong to the assigned mesh but | | 
they were assigned to the key factor 
PRESS_REP. 

To fix that we will need to get back to Eficas 
Application, locate AFFE_MODELE, delete 
GROUP_MA and assign the mesh to 
all(TOUT=OUD. 


| Ues bases sont’ fermees 


fin Ck a'execution de JOC an WXTE 


PeOUrevipeverevn rest rena rev gcan oan angen cay gegen ean gra ecea geen anne 
| eS Exception utslissteur levee sais pas interceptee. 


les 6 nailles inprinées ci-dessus n’appartiennent pas av modéle et pourtant 
| elles ont ate afvectees dane de not-cle facteur t 


Hen PeeTireag ean eoneeneenncenoenioengeennieygueanieniianiianianinanianay 


73. 


After running the Analysis, we will see that everything is going as they are 
supposed to. 

Post-Pro branch appears and your Object Browser should look like in the 
figure on your right. 


If you notice, there is one branch missing from the Fields group. 

In the IMPR_RESU command we assigned 4 different results that we want 
Code_Aster to present in Post-Pro module (DEPL, REAC_NODA, 
SIGM_NOEU and SIGM_ELNO). 

We are missing REAC_NODA, and that is why we did not ask 
CALC_CHAMP to calculate it, and IMPR_RESU to print it. 
Locate CALC_CHAMP and IMPR_RESU with the help of Efic 
Application, Find and add FORCE from the “Mots Clef Permis” list, and 
assign REAC_NODA. And do the same in RESU, IMPR_RESU concept. 


Update and re-run the Analysis. 


YY Max Elmont area 1 
Max Sie 
8 & Algonthms 
"regula 10 
1 MEASTO 20 
2 thw Mesh 1 
a eficas 
lm Aster 
Bar 
a Posteo 
Presentations 
+ Bar med 
mesh 
Families 
i) Groups 
8 Fekds 
 result_DEPL 
= sso 
i 
result SiG NGEU: 
Nodes 


© result_SIGM_ELNO, 
micelle 


i, 


The reason I am showing you all these errors is to give you a hint on how to locate your mistake, translate it, 
study it, find a solution, correct it, update the stream data on Salome_Meca and re-run the Analysis. 
Unfortunately you will spent most of your time finding and correcting your code mistakes. In this case we were 


lucky because we only made one error and it was an easy to handle, one. 


Below you can see how our command file was in the beginning, after we finished with Eficas Application, and 


how it looked after our modifications. 
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# 


Description 


15. 


eon LIRE NAILS 


soca -APFE DEL 


HAE 


(FORAT= ED" 


FCOROUE "L036"... 


soter :u0h°FE MATEREAIQULILAGEaeh 


acs FPF 


= ‘in BATERORer 30, 


result=CALC HAN 


rare Seb 


TenenaterLa,),) 


(200 "Load, 


ESL ORY 9.92 


doa 


£0, S10 20" 


1 eu 


[Figure 


#0 ere ees #2 


prcnoue nccaNzau 


rater sou-AFFE ATGRIAU( WALLACE“ 


FLAeAFFE AR MECALWODELEO Me 
otro GaOUF NOR, 


Nes HEP 


StF 200"). 
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# Enable Post-Pro module and start interrogating the results of our Finite Element Analysis. 


# |Description Figure 
1. | To acquire the deformation visualization, do the same steps as usual 
ght-click on 0,- > Deformed Shape ~> Tick Magnitude Coloring and leave the rest default 
ght-click on Def.Shape > Show Only > Representation > Surfaceframe, 

to the figure below. 


And you should have something similz 


Min. Deformation = 0 mm, Max. Deformation = 0.040 mm 


2; 
3 
result_SIGM_NOEU 
ght-click on 0,- > Scalar Map > Scalar Mode = Modulus. (default) 
ght-click on Scalar Map > Show Only > Representation > Surfaceframe. 
And you should have something similar to the figures below. 
4 BZ 
45 


6. | SIGM_NOEU ® Averaged Stress. 
Min. Stre: ).08 MPa, Max. Stress = 78.3 MPa. 
SIGM_ELNO > Unaveraged Suess. 
Min. Su LL MPa, Max. Stress 


10.7 MPa. 
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camer 


Modifying the Geometry 


# —Atthis point we are going to modify our geometry. We will use fillet on the sharp edges and create a radius, In 
order to do that, we will need to create a whole new study. Feel free to create a new folder in whatever location 
you want, and name it Filleted_Ba 

# We will use the same geometry, bar, and we are going to modify 
them into your newly created folder. 


. copy the .hdf and the comm file and paste 


Step |: Model Geometry 


# |Description Figure 

7. | Open Salome_Meca and load your -hdf file containing the previous task (Bar). 

In order to modify our geometry and then follow the steps for the FE Analysis, it would be easier for you if you 
delete all parameters below Geometry, such us Mesh, Aster Study ¢: 


2,_| Start-up Geometry Module. > 1DFillee Construction 5s 

File 10: 

Go to Operations > Fillet_ID. = 

In the pop-up window that appears, select Sketch_1 for “Wire | 

with Planar Edges”. Name [Fler 39-——————— 

Select thi th the two lines (Sharp Corner). ae 

elect the vertex that connects the two lines (Sharp Comer), 

peat a Wire with Manan Edges_@ | [= 
haere [fe [sect ietoe 

Input 15mm for “Radius”, sone 4 


F Fuse colinesr edges to allow bigger radius 


( Advanced options 
F Preview 


(aorem]_ wow |_ oe |_| 


Click “Apply and Close" 
3._ | Now go to Operations > Transformations > Mirror Image. 


vino 
In the pop-up window that appears, select the 2™ Algorithm for || ¢ 8 ef a 
“Mirror”. Sars 
: ame [oir 
Leave default name. 
St SS 
Choose Fillet ID_1 for “Objects”. oye _@ | 
ase minor [@] [Or 
Choose OX Plane for “Axe Mirror”. pas 
savanced options 
TT Set presentation parameters and sub-shapes from arguments 
1 Add pref to names of etared sb shapes 
7 Prevow 
Finish with “Apply and Close” ||| #sedcese| “aoe | cose | us| 
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Description 


Figure 


Time to create the faces out of the two drawings. 
Go to New Entity > Build > Face. 


For Face_3 choose Fillet 1D_1 for “Objects” and for Face_4 
select Mirror_2 for “Objects”. 


Finish with “Apply and Close” 


: create Face x 
Face: 


cu 


Result name 
namefces SSS 
(Face creation from wires and/or edges: 


‘Objects jefe 


7 Ty ta create 2 planar face 


(aac wy |_| aw | 


Next part is to create a shape out of the two faces we assigned. 
Go to Operations > Boolean > Fuse as we did before. 

Leave default name, and assign the Arguments accordingly. 
Face_3 for “Object 1”, 


Face_4 for “Object 2”. 


Finish with “Apply and Close” 


= 
ea 
————$<$< 
SS 
segues 
‘Object 1 el es °° 
oct [@ [x7 


danced sptions 


TT Set presentation parameters and sub-shapes from arguments 
IF Add profi to names of restored sub shapes 


Now you have created the two fillets needed for our Analysis. 


Save your work, be careful not to overwrite with the previous 
task, and let’s continue to assign our groups. 


You don’t need to delete the old groups which we assigned for 
Fuse_l. 


We will create the same groups for Fuse_2 now. 


Your Object Browser should look like the figure on the right. 


(lwescme) wow | coe | _ ute | 
z 


‘Object Browser 


= Sketch 
Mirror_2 
Face 3 
Face 4 
Hiee's 
Face 3 
«Face 4 
BZ) Fx 

BZ) Load 
es 


cea 
5 


‘Guide 


Save again, and open Mesh Module from the drop down list. 
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Step I: Meshing Geometry 
# |Description Figure a 
7. | As soon as Mesh Module loads, go to Mesh > Create Mesh, - Ee 

and fill in the blanks like the picture on the right. cers eh? 
Geomet Fuse? 
The Algorithms and the Hypothesis for both 2D and 1D tabs Hs 2 
are the same like in the non-filleted geometry. 
y so 2] w| oo] 
Furthermore, Salome_Meca has stored the Hypothesis Algorithm — iangie (meio) Sd 
Parameters so you won’t have to assign them again. oe. TZ | 
30 | 2 10] o] 
Agorthm [Wire Decretzaten >] 
Finish with “Apply and Close” || Hypothesis [WaxSeea SSCS 
2._| Right-click on Mesh_2 > Create Groups from Geometry. 
oer ee 
Right-click on Mesh_2 > Convert to/from quadratic. ca 
ame 
Right-click on Mesh_2 > Compute. wash 2 
shes 
In the figure on the right, there is my mesh information. Ta! thnear uedratc 
Nodes: en 
We did not do anything different than the previous Mesh_1. pene 
alle ° 
ges: " ° e 
recast os ° sa 
ingles 534 0 4 
oustranaies! a ° ° 
Poygons: 0 
whines: oo ° ° 
ishedrone: 8 0 ° 
Hevahesrone: 0 ° ° 
Pyramids ‘ ° ° 
Poms: ° ° ° 
Menageal prams: 0 
Poiyredrne= 0 
3. | That would be all the necessary information to modify for both Geometry and Mesh Module. Let’s continue now | 
to create our command file. In this case we will not have to change anything at all to our command file, because 
the changes we did in the Geometry are not affecting the commands. Also we have the same groups from 
Geometry who also got the same names. 
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#_| Description 


1. | Open up Aster Module from the drop down list. 
Go to Aster > Add Study Case, 


Fill in the gaps according to the figure on the 
right. 


Notice that from Mesh, I chose Mesh_2. Which 
means that even though we got the same group 
names in Fuse_| and Fuse_2, the Analysis will 
take into consideration only the geometry which 
we assigned to the Mesh_2, that would be 


ho 


Finish with “Ok” 


[Figure 
aesubppiation 
sy cave date 
mz 
Cammand ie [For dek =] ah frees. sarracenarcar comm 
Mesh [fom spect browser Ie _|fmena 


AST serves: 


Server [[ecamos x] Aserveron [STAI =] _lesh servers 


Execution mode nTe=eT —T] Interactive fotow up FF 


Suber prametrs 
‘etal memory (me) [ics] __ times) [zo 
(CU number [7 Sy save rma database 


2._ | Now itis time to run the Analysis. 
Right-click to “Firstry” (or the name you ass 


ned for the study), select Run. Everything should work out fine, 
without facing any problems or errors, and Post-Pro branch should appear. 


CHAP! 
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# |Description Figure 
1. | To acquire the deformation visualization, do the same steps as usual 
t-click on 0,- > Deformed Shape > Tick Magnitude Coloring and leave the rest default 
ght-click on Def.Shape > Show Only > Representation > Surfaceframe, 
And you should have something similar to the figure below 
Min. Deformation = 5.56e-22 mm, Max. Deformation = 0.039 1 
2. 
3. | To check the Normal Avera; iged Stress results, we have to open both result_SIGM_ELNO and 
result_SIGM_NOEU 
{-click on 0,- > Scalar Map > Scalar Mode = Modulus, (default) 
t-click on Scalar Map > Show Only > Representation > Surfaceframe. 
And you should have something similar to the figures below. 
4 s 
5 
6 
Min. Stress = 0.28 MPa, Max. Stress = 74.96 MPa. 
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Finite Element Convergence 


Plate Bar: 
One of the core concepts of the finite element method is that, as mentioned, the finer the mesh, the more accurate 


the solution, Ultimately, the solution will reach the analytical solution. But how fast does it approach that solution? This is 
what we want to answer in this part of the chapter. 

The answer depends on what kind of element we are using. In this step we will add a Sub-Mesh to our model, to 
make the mesh finer in a certain point of interest, and try to compare different kinds of Sub-Mesh. The comparison will of 
course be in the same part of the Geometry so that the results are comparable. 

# To do so, we must go back to the original Plate Bar and make some modifications, after that we will modify also 
the Filleted Bar. 


#_| Description Figure 
T._| Open Salome_Meca and load the Plate Bar case study. : 
Sree 
Activate Geometry Module, in order to create some more ae Go a ce 
groups to help us create the Sub-Mesh. I ——— eee 
Name [etek 
Select Fuse_l, and go to New Entity > Group > Create. 
‘ain Shape nd Subabapes 
We will make two different groups, one will contain the two Pes sis | ve 
vertical lines, and the other one will contain the horizontal ‘ein Shope Sletion retrision show only sete 
lines connected to those vertical ones. Laie sive) Hide soot 
(© Geometries pars eft second shape | EES _I 
© oniysusshapes the Second shape | SOW i sb-shope| 
second Stace _@ | 
=—__—_—— Select All 
1 
Shove _| _ Help 


2._| After creating the groups, activate Mesh Module from the drop Crete se 7 
down list. [sume 
Go to Mesh > Create Sub-mesh, and fill the gaps according to a 
the picture on the right. I have named VertMesh_1 and ‘Geometry [ @ [veroresh 2 

HorMesh_l the four Vertical and Horizontal lines we are 


sp] 20 19} oo] 


interested in making a finer mesh. 
iypothess = alal 
‘Add. Hypothesis [=tione> a 
‘Asigi a sot of iypotiosse | 
‘Avoly and close] ply | chose |__eip 


3._| This figure shows the Hypothesis Argument/settings. SSS 
The aim of this section is to calculate different results by ‘umber of segments [3a] 
messing with the mesh. Pees [sae dntn il 
More specific, we will divide our Vertical lines to 5 (Number | sealermtoe [OSS 
of Segments) elements. 

But we will use Scale Factor of 0.05 to distribute the elements 
differently between the lines length, making them finer and 
smaller as we reach the connection point. 
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Hypo 


#_| Description [Figure 
4._| For the Horizontal lines, my Hypothesis 
Arguments are similar to the figure on the Name Fey ~ ‘iypothessConetiction 
right. nico | Fan 8 Number or segments 
Number of Segments = 10, with Seale Factor | | somsey Pw [fau==ar | “Une 
of 5. Number oF Segments [10 
sp} 20 0 | oo] || Woeetabuinwon [Scaeacumuto 
thm [Eempta | Stal Factor a 


‘Add. Hypothesis 


‘Apply and cosa 


Assia 


| _ cancel 


# We know now how to create a Sub-mesh. I will summarize one more time, Sub-mesh gives us the ability to point 
out a specific area in our geometry and input different Mesh options like the ones we got in the original Mesh. 


# Finite Element Convergence is all about figuring the right Mesh and Sub-mesh, to get results as closer to reality 


as possible. 


# To achieve that, Mesh needs to get finer. So we will adjust the Number of Segments and the Scale Factor of the 


Sub-Meshes in order to analyze the results, 


Results 


Mesh: 2D 
Algorithm 
Hypothesis 


Tam using > 


[VertMesh] SubMesh_I: 1D 
Algorithm 


Hypothesis 


> Triangle (Mefisto) 
> Max. Element Area_| 


Arguments 


[HorMesh]  SubMesh_I: 1D 
Algorithm 


Hypothesis 


Arguments 


1D 
Algorithm > Wire Discretisation 
20. Hypothesis. > Max. Size_| = 10 


> Composite Side Discretisation 
> Nb. Segments_1 


> Number of Segments = 15 
Scale Distribution 
Seale Factor = 


> Composite Side Diseretisation 
> Nb. Segments_2 


> Number of Segments = 25 
Seale Distribution 
Scale Factor = 20 
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# These are my starting configurations for the Mesh, Sub-meshes. 


# After you have input all the Algorithms/Hypothesis/Arguments, you must now get started to calculate the given 
data. 

# The procedure that you will follow to manually change the Arguments is actually very simple, but time 
consuming. 


1, Activate Mesh Module. 

2, Right-click one SubMesh_1 and choose Edit Mesh/Submesh. 

3. Get to Hypothesis Construction, when the pop-up window appears. 
4, Change the values and finish with “Apply and Close” 

5. Do the same for SubMesh_2. 

6. Right-click on Mesh_1 and choose Compute. 

7. Right-click on Mesh_1 and choose Conyert to/from quadratic. 

8. Right-click on Mesh_1 and choose Compute. 

9. Activate Aster Module. 

10. Right-click on “assignedname” study and choose Update Mesh. 

11. Right-click on “assignedname” study and choose Run. 

12. Activate Post-Pro Module after the Analysis is done. 

13. Get to the Averaged Stress tab, that would be result_SIGM_NOEU. 
14, Assign Scalar Map with default options to 0,- . 

15. Right-click on Scalar Map and choose Representation > Surfaceframe. 
16. Right-click on Scalar Map and choose Show Only. 


# Here are my own results, 


Nb.Seg/Scal, Fac. _Nb.Seg/Scal. Fac. Nb.Seg/Scal. Fac. Nb.Seg/Scal. Fac. _Nb.Seg/Scal. Fac. 


ala’ z 2 3 4 5 
Sub-mesh_1 15/1 15/0.8 15/0.6 15/0.4 15/0.2 
Sub-mesh_2 25/20 30/20 35/20 35/30 35/33 

Averaged Stress 189.628 197.76 212.436 289.679 304.92 


tarting Mesh, 2 = Good Mesh, 3 = Very Good Mesh, 4 = Excellent Mesh, 5 = Ultimate Mesh. 
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# You can see from the results of Max. Averaged Normal Stress, that from #1 to #5 (the mesh is getting finer) the 


difference between Stresses increases. 

# For example 
Second Difference [ 212.4 
Third Difference [ 289.7 -21 
Forth Difference [ 304.9 


197.8 


First Difference [ 197.8 - 189.6 = ~8 MPa] 

15 MPa] 
77 MPa] 
289.7 = 15 MPa] 


# We can understand from those results that, the finer the mesh gets, the difference between the two results is 


getting bigger, theoretically touching infinite. 


Filleted_Plate_Bar: 


# After finishing up the simple Plate Bar, it is time to do the same to the Filleted_Plate_Bar. 


#_| Description Figure 7] 
7. | Open up Salome_Meca and load the Filleted_Bar study. : creat croup - 
Shape We 
Activate Geometry Module, ct a ca 
‘Group Nore 
Create a group of the fillets this time, as shown in the figure on | || wane [Flas 
the right. ‘Male Shape Aad Sibeshapes $$ 
| 
= Man shape selection reticion Show only selected 
¢ mvc ners | 
© ceometical pare fthe second shape === _] 
© cntysubshapes ofthe secon shape | Show a sbshapes| 
second Shape @ | 
e4 Select AL 
ii 
Finish with “Apply and Close” | | |sz=eiee=a] sso _| gee | _ Ho 
2. | After creating the groups, activate Mesh Module from the drop || Edit meshsub mosh | 
down list. name [F eon s._<— T 
en f= 
Go to Mesh > Create Sub-mesh, and fill the gaps according to | |eemety [Frets | 
the picture on the right, 
30] 20 | 0] 
This time we will use Arithmetic 1D_1 for “Hypothes Agorthm [Compote Side Decreteaton =] 
Hypothesis [arthmase D1 =] asl al 
Afier assigning the Hypothesis don’t forget to input the data fd, typothess [nes SSCS 


for the Arguments, 


After that, Finish with “Apply and Close” 


3._| This figure shows the Hypothesis Argument/settings. 
Start Length is the length of our first element, 

End Length is the length of the last element. 

The length from one element to the other is increasing 
automatically, starting from 0.2 up to 3. 


‘sign ae athypotnenne J 


eguments 
stactenetn [EE 
a) 
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# Following the same procedure like we did in the Plate_Bar, you should have your results in a couple of minutes. 


Results 


Tamusing > Mesh: 2D 1D 
Algorithm > Triangle (Mefisto) Algorithm > Wire Discretisation 
Hypothesis > Max. Element Area_! = 20. Hypothesis > Max. Size_1 = 10 


[Fillets] SubMesh_l: 1D 
Algorithm 


Hypothesis 


# Here are my own results: 


Nb.Seg/Scal. Fac. 


> Composite Side Discretisation 
> Arithmetic 1D_1 
Arguments > 


Start Length 
End Length 


Nb.Seg/Scal. Fac. Nb.Seg/Scal. Fac. Nb.Seg/Scal. Fac. 


ala? 1 2 3 4 
Sub-mesh_1 0.2-3 0.2-2 l 02-1 0.2-05 
Averaged Stress 79.842 82.738 84.155 84.417 


# Here, we can observe a different phenomenon than the one we 


\w before. In this case elements of the Sub-mesh 


are getting finer throughout the fillets length the Averaged Stress result difference is getting smaller, therefore 


convergent to some precise value. 


starting Mesh, 


= Very Good Mesh, 4 = Excellent Mesh, 5 = Ultimate Mesh. 
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Chapter_VI: Shaft 


In this Chapter we are interested in a 2D Axisymmetric 


Analysis with a notched 


2D_ANALYsIS, AXI_SYMMETRY, NOTCH 


Step 1: Purpose of the FE Analysis/ Description of the problem 


#= = Axisymmetric Analysis - The axisymmetric problem deals with the analysis of structures of revolution under 
axisymmetric loading. A structure of revolution is generated by a generating cross section that rotates 360° about 
an axis of revolution, Such structures are said to be rotationally symmetric, 

# A structure of revolution by itself does not necessarily define an axisymmetric problem. It is also necessary that 
the loading, as well as the support boundary conditions, be rotationally symmetric. 

# Axisymmetric elements are 2D elements that can be used to model axisymmetric geometries with axisymmetric 
load. In simpler words, we are converting a 3D Geometry to a 2D Geometry making the model smaller, therefore 
faster execution and faster post processing. 

# We only model the cross section, and Salome_Meca accounts for the fact that it is really a 3D, axisymmetric 
structure. 

# The axisymmetric problem deals with the analysis of structures of revolution under axisymmetric loading. A 
structure of revolution is generated by a generating cross section that rotates 360° about an axis of revolution, 

# In this Chapter we consider the finite element discretization of axisymmetric solids. We were given a Steel Shaft 
with dimensions which are given below. Purpose of this chapter is to fully understand the concept of Axi- 
symmetry and its possibilities, 


Examples before starting ourmodel. = === 


fy | € Surfaces Needed. CoD 
Models > 
# In Salome_Meca, 
axisymmetric models must 
be drawn in OY Plane. 
# The OX direction is the 
radial direction. 
# The 2D model will be rotated = eS 
about the Y axis (always i 
L\ 
| 


about x=0). 


# Nothing in your model 
should be in the region X<0, 


# In post processing o, will be 
the radial stress, oy will be 
the axial stress, and 0; will be 


_ lana 
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Step 2: Input Values for the FE Analysis 


Shaft Dimensions: 


Young’s Modulus: 
Poisson’s Rati 
Pressure Applied: 


Length 


Length = 75 mm [1]; 
Radius = 5 mm [R]; 
200 GPa; 

03; 


50 MPa [P]; 


50 mm [h]; 


Length = 15 mm [d]; 
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Step 3: Model Geometry 


2D_ANALYsIS, AXI_SYMMETRY, NOTCH 


# Start Salome_Meca, create a new document, save as “Shaft”, activate Geometry Module from the drop down list 
and start sketching. 
# For generating an Axi-symmetric model, we just need to create a planar face of the meridian cross section of the 
shaft. 
#_| Description Figure 
T. [In the Geometry Module, go to New Entity > Basics > 2D . 20 Skatch Construction x 
Sketch. Coordinate system 
[abet coovinate atom =] | Restore 
These are the points we are going to use to create our 2D 
gece —k——— 
: "Element ype 
Point_I: (15,0) Point_4: (20,75) ae ees ir 
Point_2: (15,150) Point_S: (25,75) 
Point_3: (20,150) Point_6: (25,0) coke 
5 —— 
© Point © Direction 
Fone Aatitonnal Parameters 
© Absolte Nona Tangent) 
© Rehtve © fone 
© Selection © center 
values 
xfs = eo 
mg A [ono 
Finish with “Sketch Closure” Close _| Sketch Gosure| cancel | Hen 
2._| After creating our 2D sketch, it is time to create the radius. : 1p Fillet Construction 
Fille 10- 
Go to Operations > Fillet 1D. ey 
Select with the mouse the correct Vertex and input Smm for esate nm} $$ > 
Radius, according to the figure on the right. Name [Filetip. 
file 1D 
Wire with Panar Edges _@ | [Sketch 1 
‘Vertexes @ | [stetch aventexe 
ads: | 
IF Fuse collinear ages to allow bigger radius 
Advanced options 
F Proview 


Auply and Clse| Apply | Close Help 


Finish with “Apply and Close” 


ep cuaprer_V! 


HAF 
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# |Description Figure 1 
3. | Next step is creating a face out of the sketch, in which we will || create A Face x 
assign the groups later on, Face 
re 
Go to New Entity > Build > Face. 
es pam 
Select Fillet 1D_I for Object. mame [Faced 
Face creation fram wires and/or edges 
‘Objects | @ | [Fer 1.2 
F ny tocreate a plana face 
toy | one |e | 
Finish with “Apply and Close" 
4. | After creating a face, we are going to assign the groups that we | = creat Group ra 
are going to need later on, Shape pe 
rie 6 a 
Go to New Entity > Group > Create. roup wae 
a 
Select 2 Algorithm for Lines, name the first group “Fix”, and 
s a 
choose for Main Shape “Face_1”. 
Select the Line where we want our Fixed Support to be, and Main shape _@ | [Face 
press “Add”, Main shape selection etcon Sr 
© Noresiration er 
Do the same for the load. Create a group named “Load” with cS Sas Shaner eas 
the same Line Algorithm and “Add” the Line where our loads 3 lt | aera 
ie Secondshape @|T 
sett 
Also you can create one more group for the fillet, if you want — 
to make the mesh in that part finer. = 
Now for the “Guide” group, select 3“ Algorithm for fe 
“Add” the only face available to you. opty nd Cove] _toply Ser | 
Finish with “Apply and Close” 
5. | Off-Topic. > Constrition Sy Revelation x 
Revolution 
Another cool ability of Salome_Meca is the Revolution option. ||| ¢ gy 
In this task we were asked to create a face which if you revolve | | ; Bewtnane 
it 360° it would turn into a shaft. ame [Revelan 
That is what the Revolution command does. pes 
objects _@ | [Fe 
Go to New Entity > Generation > Revolution ce |p 
fi, a SF Reverse’ 
Input the data like the figure on the right. seis — 
T sath Directions 
(Advanced options 
F rrevew 
Finish with “Apply and Close": 
ee and ios Jose us 
‘And you can see now in the Viewer tab, your half-Shatt. seers | apo | cere |e] 
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Step 4: Meshing of Geometry 
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# After saving your geometry study, activate Mesh Module from the drop down list and start meshing your model. 
# |Description Figure J 
1. | Highlight Face_I, and go to Mesh > Create Mesh. : create mesh +e 
ame [re 
Leave default name “Mesh_1” and fill in the blanks according | semen [@ [<== 
to the figure on the right. 
20 20] 1} 0] 
For Hypothesis parameters I applies for Max. Element Size= 3, | Asati — regent02 ad 
and Min, Element Size= 2. Viypotness [RETGENDOTeane=s __=] ai) BI 
‘dd. ypothess [ete a 
“assign a set of hypotheses J 
Finish with “Apply and Close” | eeijasa isi soni cose tip 
2._| Next part is creating a Sub-Mesh in the filleted area. + Create sabmesh om 
name [wen 
Go to Mesh > Create Sub-Mesh, Meh |= 
eometry [@ [Fe 
Fill in the gaps according to the figure on the right, BUT be 
careful for Geometry input you have to choose the group we =| 20 | op] 
made before that includes the fillet edge. Agaiim — [Wrebaramaion =] 
Hypathess = Masel Sa 
Max. Element Size = 0.5 a 
sigh aetna | 
Finish with “Apply and Clos sooty | chose |] tie 
3. | After computing the Mesh_I which includes SubMesh_I, right [/~ Nosh computation suceed 7 
click on Mesh_I and choose “Convert to/from Quadratic mpute mesh $$} 
+e 
Right-click on Mesh_1 again and choose compute. Name 
Mesh 2 
The figure on the right shows my Mesh Information. Mech nior 
al ner ude 
Nodes an 
‘OD Eloments: 0 
ale: ° 
ages @ aa 
Facer? @ we 
‘iengles 08 ° 08 
uacranges: 0 ° ° 
foygene: 8 
Volumes: fy ° 
marhedrns: 0 ° ° 
Herohetrons:— @ ° ° 
ryamiss: 0 ° ° 
Hexaponl pss: 0 
Poypedions: 0 
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Step 5, 6, 7, 8, 9, 10: Salome_Meca Linear Static Analysis Wizard, Running the Analysis 


# After Meshing your Geometry, it’s time to create the command file. Activate Aster Module from the drop down 
list and click on Aster > Wizards > Linear Elastic. 


#_| Description Figure 
7._| In the model definition window that pops up select “Axis symmetric” and click “Next”. 


Select “Mesh_1” in the next step if not already selected and click “Next”. 


Enter 2e5 for Young’s Modulus and 0.3 for Poisson’s Ratio and click “Next”. 
2._| For the Geometric Boundary Conditions, select “Fi 


for the Group tab. teotropic linear elaate study 
undies condone 


: ‘ot-subapplication = 


Delete “0” from “DX” as we do not want to constrain 


the Shaft in X direction. _& | Adding imposed degrees of resdom on groups 
For “DY” leave it default as we want it to not move in Sroup on Ls 
Y direction. BI L 
tek cancel 
Proceed with “Next” 
3._| For the Load Boundary Conditions, select “Load” for |~> eaubappiaton x 
the Group tab. Isotropic linear elastic study 
oundaves conditions 


Enter *-50" as “Pressure” 


4 | ing pearson mee sop 
Proceed with 


Group Pressure 


[se 


Enter a file name, and save the command file to the 
same folder you saved your geometry study. 


Finish with “Finish” option in the wizard window. 


<Back 
4_| Check that a new branch of “Aster” has been added in the Object Browser. This concludes setting up the study 
for Finite Element Analysis with Salome_Meca Wizard. 


Right-click on “linear-static” branch in Aster and select “Run”. 


If everything goes well a new branch of “Post-Pro” will be added to the Object Browser. 
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it’s time to analyze the results. Activate Post-Pro Module 


After running the analysis and Post-Pro branch appe: 
from the drop down list and 


Description Figure 


“Deformed Shape” and 


Open branch “RESU__DEPL” and right click on From the pop up window sele 
visualize it 


w how much it has deformed, right-click on “Mesh_1” in the Mesh Module without activating it and select 


Do the same for branch “RESU__SIGM_NOEU” and“RESU__SIEQ_NO) 
For the Normal stress (SIGM) in the pop up window change Scalar Mode to “SIY 


symmetry) 
For the Equivalent stress (SIEQ) in the pop up window chang. 


*. (Y axis direction for axi- 


Scalar Mode to “VMIS” 


Thave placed my own results in the figures below. 


RESU____DEPL 0, - ‘RESU___SIGM_NOEU 0, - RESU___SIEQ NOE 0, - 


0.0271814 | 84.3782 77.0214 


0.0138368 41.6371 


20.2665 23.5428 


0.000492144 -1.10401 
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Stress Concentration Factor 


# A stress concentration is a location in an object where stress is concentrated. An object is strongest when force is 
evenly distributed over its area, so a reduction in area, e.g. caused by a crack, results in a localized increase in 
stress, A material fail, via a propagating crack, when a concentrated stress exceeds the material’s theoretical 
cohesive strength. The real fracture strength of a material is always lower than the theoretical value because most 
materials contain small cracks or contaminants that concentrate stress. 


# The stress concentrators are geometrical irregularities that cause an increase in the average effort that should be 
present in regions near these discontinuities, the relationship between the maximum stress that occurs and the 
average effort that should occur is defined as stress concentration factor; which is determined by experimental or 
analytical methods and presented in graphical form for ease interpretation, 


# The stress concentration factor for a tube in tension with fillet, our case, can be determined as the relation of the 
maximum normal stress in the discontinuity and the nominal stress, and is obtained through the equation: 


ToT 
tlh = 0.25 
2.0 
}~ 0.50, 
18 / 4,00 
1.6 3.00, 
RK, 
1.4) 
1.2) 
1.0' 
0 1 2 3 4 
tir 
Chart Stress concentration factors K, for a tube in tension with fillet (Lee and Ades 1956; 


ESDU 1981). 


eyyycuaprer_v! 
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Hand Calculations VS Computational Calculations of Stress Concentration 


Hand Calculations 


As we can see in the graph, there are 2 parameters that are taken into consideration to aequi 


re stress concentration 


factor. First parameter would be to determine which line do we need to choose, to accomplish that we have to solve the 


equation t/A with our given values. Second parameter would be to acquire the correct output for the X axi 


equation t/r . 


According to the graph, my results are: . 


Computational Calculations 


olving the 


. Ki= 1.52) 


To acquire the Stress Concentration Factor for the computational calculations we are going to need Max. Normal 
Stress (SIGM_NOEU > SIYY) and Nominal Stress, which in our case equals to Pressure = SOMPa. 


We don't need to calculate the Max. Normal St 


ince Salome_Meca has an option of showing it to us, 


#_| Description 


Figure 


T._| This is a good opportunity to check another aspect of 
Salome_Meca’s Post-Pro Module, an easy way to pick points 
throughout a geometry. 


Activate Post-Pro Module from the drop down list, open 
RESU__ SIGM_NOEU branch and highlight “Scalar Map: 


In the top right corner, there is a group of 5 icons, identical 
with the ones shown in the figure on the right. 
Activate “Selection Panel 


Bvls & 


SELECTION PANEL 
ACTOR SELECTION 
CELL SELECTION 
POINT SELECTION 
GAUSS POINT SELECTION 


2._| When you activate “Selection Panel”, it will open up a new 
window on the right side of Salome Platform. 


That window has 4 main tabs. 
[Actor, Cell, Point, Gauss Point]. 


The tab that we are interested in this, 
case is the “Point” tab. 

By choosing “Point” tab and by 
pointing out that we are interested 
for the point where the Scalar Value 
s “Maximum”, the software will 
give you all the information 
regarding the Maximum Scalar 
Value Point which in our case is also 
the Maximum Normal Stress Point. 


4 
SSS 
os 
asi 
== f 
es 


3. | According to the computational results: Ky 


84.3782.) => K,~ 1.69 


Solving the Equation 


© (di/h + d;/t) >28 => (30/5 + 30/5) >28 =>12>28 @ 


Conclusion: As we can see the results between Hand and Computational Calculation 


are almost the same. The 


difference of 11% is probably because of the geometry. Geometry proportion does not meet the criterion for the 
experimental data in the Stress Concentration Factor diagram, [Eq. @)]. The Geometry proportion coefficient should be 


ater than 28, 


g 
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Redefining Mesh 


# We can see from the Post-Pro Module information results, that our mesh is a little irregular at some points, that is 
a very good reason to refine Mesh to almost half of its current size, and see what happens to our results 


# [Description Figure 
T. | In order to do so, we need to activate Mesh Module again. Tp ietiels conerion oe 
1 wager 20 
After activating Mesh Module, open Mesh branch, open Hypotheses branch || arguments | tscal sees | 
and right click to NETGEN 2D Parameters and reduce Max. and Min. Size to || yy soe ——s 
half of the current value. rae = 
Peseeores 
Do the same for Max Size_1 Hypothesis, and reduce to half the Length value, || [=m 


Right click on Mesh and Re-Compute it. oes aa =| 

Right click on Mesh and Convert to/from quadratic. pesreas |B EI 

Right click on Mesh and Re-Compute it. Nb. sogspernadus [E—__] 
TF alow quasrangies 


1 Optimize 


| _ cancel He 


2. | Next step is activating the Aster Module, in order to re-calculate the new adjustments. 


Right click to “linear-static” and choose Update Mesh. 
Right click to “linear-static” and choose Run. 


3._| Post-Pro branch is now visible. Activate Post-Pro Module from the drop down list. 
Open the branches until u reach RESU__SIGM_NOEU’s “0,-* and right click on it. 
Choose Scalar Map with Scalar Mode = SIYY. 


Right click on it and choose Show Only. 
4_| The figure on the right, indicates my own results. 


84,3428 


We can see that the difference between these results and the previous ones is really small, but it is 
obvious that when we you change Mesh Parameters the results are changing also. 

62.9930 
In this case the difference is a decrease of 0.05%. 

A decrease like that, can hardly change the Stress Concentration Factor. 


K, = 84-3428/,) => K,~1.69 


1.05644] 


CHAPTER_VII: LEVEL_OF_GEOMETRY 
Chapter_VI: Level of Geometry 
In this Chapter we are interested in comparing 1D, 2D, 3D 


Analysis methods. 
Beams, Solids and Shell Element Method. 


aod 


1D - 2D - 3D_ANALYSIS_COMPARISON 
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Step 1: Purpose of the FE Analysis/ Description of the Problem a 


# Beam Element - A beam element is a slender structural member that offers resistance to forces and bending 
under applied loads, A beam element differs from a truss element in that a beam resists moments (twisting and 
bending) at the connections. 


# These three node elements are formulated in three-dimensional space. The first two nodes (I-node and J-node) are 
specified by the element geometry. The third node (K-node) is used to orient each beam element in 3-D space 
(see Figure 1). A maximum of three translational degrees-of-freedom and three rotational degrees-of-freedom are 
defined for beam elements (see Figure 2). Three orthogonal forces (one axial and two shear) and three orthogonal 
moments (one torsion and two bending) are calculated at each end of each element. Optionally, the maximum 
normal stresses produced by combined axial and bending loads are calculated. Uniform inertia loads in three 
directions, fixed-end forces, and intermediate loads are the basic element based loadings. 


# The basic guidelines for when to use a beam element are: 
© The length of the element is much greater than the width or depth, 
© The element has constant cross-sectional properties. 
© The element must be able to transfer moments. 
© The element must be able to handle a load distributed across its length, 


# Solid Element ~ Solid elements are three-dimensional finite elements that can 
model solid bodies and structures without any a priori geometric simplification. 


# Finite element models of this type have the advantage of directness. Geometric, 
constitutive and loading assumptions required to effect dimensionality reduction, 
for example to planar or axisymmetric behavior, are avoided. Boundary 
conditions on both forces and displacements can be more realistically treated. 
Another attractive feature is that the finite element mesh visually looks like the 
physical system. 


# Summarizing, use of solid elements should be restricted to problem and analysis gfe 2:Boam Elomont Dogrove ot Freedom 
stages, such as verification, where the generality and flexibility of full 3D models 
is warranted. They should be avoided during design stages. Furthermore they should also be avoided in thin-wall 
structures such as aerospace shells, since solid elements tend to perform poorly because of locking problems. 


# Shell Element — Shell elements are 4-8 node isoparametric quadrilaterals or 3-6 node y " 


triangular elements in any 3D orientation. The 4 node elements require a much finer 4 ee 
mesh than the 8 node elements to give convergent displacements and stresses in ‘a 
models involving out-of-plane bending. Figure | shows some typical shell elements. 

\w 


# The general and co-rotational shell element is formulated based on works by Ahmad, 
Iron and Zienkiewiez and later refined by Bathe and Balourchi. It can be applied to 
mode! both thick and thin shell problems. Also, the geometry of a doubly curved shell 7 

ariable thickness can be accurately described using the shell element. , 


# The thin shell element is based on thin plate theory. The bending behavior of the x 
element is based on a discrete Kirchhoff approach to plate bending using Batoz’s A LAr 
interpolation functions. This formulation satisfies Kirchhoff constraints along the 
boundary and provides linear variation of curvature through the element. The a ae 
membrane behavior of the element is based on the Allman triangle which is derived ant ant 
from the Linear Strain Triangular (LST) element. A general curved surface is 7 {4 
approximated by this element as a set of facets formed by the planes defined by the , 
three nodes of each element. For these reason a well-refined mesh is necessary x | 


# The element geometry is described by the nodal point coordinates. Each shell element 
node has 5 degrees of freedom (DOF), three translations and two rotations. The 
translational DOF are in the global Cartesian coordinate system. The rotations are about two orthogonal axes on 
the shell surface defined at each node. The rotational boundary condition restraints and applied moments also 
refer to this nodal rotational system. The two rotational axes (V1 and V2) are usually automatically determined by 
the processor and you do not have to specifically orient them, 
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Step 2: Input Values for the FE Analysis 


Chassis Dimensions: — Length = 115 mm [h}; Length = 12.5 mm [a]; 
Length = 100 mm [I]; Length = 1200 mm [m]; 
Length = 20 mm [d]; 

Young’s Modul 200 GPa: 

Poisson’s Rati 03; 

Pressure Applied: 10 MPa [P]; In case 


Supports: Fixed Support at “A” and “B” point. 
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i. Beam Elements 


Step 3: Model Geometry 


1D - 2D - 3D_ANALYSIS_COMPARISON 


# Open up Salome_Meca, activate Geometry Module from the drop down list, save your study to a desired location 
and get ready to start modeling. 
# |Description Figure | 
1. | Creating our Beam should be an easy task. ~ Point Construction =] 
Paints 
Go to New Entity > Basic > Point. es. re -O cx Co 
Create two points in order to generate the line later. ieee 
name [F2_s00,00 
Point_I = (0,0,0) Coordinates 
Point_2 = (600,0,0) *Symmetric Model* x: [500 
ye 
Finish with “Apply and Close” 2:6 
soy |__| 
2. | Next step, is creating our Line. | 
tine 
Go to New Entity > Basic > Line. ona ce 
——————— 
Insert P1_0,0,0 for the starting point. ane [== 
P2_600,0,0 for the end point. = 
——————— 
Finish with “Apply and Clo: Pune _@ | Faa8 
pont 2 [@ [F2600.0.0 
Maply end cese| apply | cose | wep | 
3._| Final step to create our Geometry, would be adjusting our = create rowp x 
groups. ee 
oe e a ca 
Go to New Entity > Groups > Create. p60 ame 
ame [Fx 
In the pop up window, select 1 Algorithm for Shape Type of | | -wainshapeand subshapa 
Nodes. Main Shape xa] fivex 
~ Nan shape Seaton reson aaay aaa 
Enter Fix as a name, for the Fixed Support point, and “Add” © io restrlction Saale | 
the Point_l. eercaremcnee—s [EE 
Click “Apply” ||| | Conrsubstaresertne sondage | Storateustans| 
Do the same for the other end point, where we want our ‘second shape _@ | 
Displacement Support to be, and name it “Disp”. _—_—__ 
Click “Apply” — 
Lastly, select the 2" Algorithm for Shape Type of Lines, name = 
the group “Load” and “Add” the only line we got available. 
Finish with “Apply and Close” 
inish with “Apply and Close” | [aes | ernen] ee ee 
4_| These steps conclude out Geometry creation session. Save your project and continue to Mesh creation. 
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Step 4: Meshing Geometry 
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# |Description Figure 
7. | Activate Mesh Module from the drop down list. ss cratemeh = 
ame [re 
Go to Mesh > Create Mesh. Geometry [@ [Une 
Select Line_I for Geometry. 30)| 205 15)] 00] 
Select Wire Discretisation for Algorithm. Agorthm [Wee Decreteaton =] 
Select Nb. Segments_1 for Hypothesis, with number of Hypothesis [NE Seqrete 1 =] a) | 
Segments “9” (it might seem strange, but we will explain why || i49 ypenese===-———] al al 
we chose 9). 
Asign 2 set of hypotheses | 
Finish with “Apply and Close” ‘Awply and Cose| apply | close |Help 
2._| Right-click on Mesh_I > Compute. sh computation sczend = 
compute mash 
If everything is done correctly, you will have the same results 
as the figure on the right. 
“Hypothesis: Nb, Segments > The ability of this hypothesis is 
to separate your geometry to as many parts (in our case Edges) Linear vada 
as you define. 
We set the number of segments to 9, which means that our 
mesh is going to consist of 9 Edges and 10 Nodes. 
> ° 
° ° 
° ° 


Click “Close” 
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Step 5, 6, 7, 8, 9, 10: Creating the command file, run the analysis 


# In order to create our command file, this time we will use Eficas MultiCatalogue, which can prevent us from 
making grammatical mistakes in Code_Aster. 
# |Description Figure 
1. | Goto Eficas > Eficas MultiCatalogue. 
In the pop-up window, click on Create New or Nouveau, 
Choose ASTER > “Ok”, STAI] > “Ok” and you have set up Eficas and it is ready for modification. 
Before starting editing, save the command file in the same location where the geometry/mesb/.hdf files are 
currently located, with an appropriate name. 
Fishier > Enregistrer sous. 
We will explain step by step the commands we are going to use to complete our command file, because it is 
going to be a little complicated, | 
2, | Beam.comm > Command file name, 
DEBUT > Start. [concaptvaleur 
DEFI_MATERIAU > Define Material. 
ELAS > Elastic (Constitutive law). aes 
E> Young’s Modulus. 
NU > Poisson’s Ratio. eam 
LIRE_MAILLAGE > Reading the mesh. mesh 
MAILLAGE > Mesh. MED 
CREA_GROUP_NO > Create Nodal group. mesh 
GROUP_MA > Group of elements. mesh 
ae 
x (Lose) 
@ b nom group_ma 
3. | AFFE_MODELE > Assign to the model. Beam.comm | 
PHENOMENE ~ Phenomenon. ennai ERC 
MODELISATION > Modelization. = 8 gamcenn 
si fetia ym DEBUT 
AFFE_MATERIAU > Assign the material. Pim inee oicianils ow 
TOUT > All. URE MAILLAGE = mesh 
DEFL_GROUP mesh 
model 
mesh 
S@ aFFE 
© Grour va (woad) 
© PHENOMENE : MECANIQUE 
© @ b meconique 
MODELISATION : (ou 
(6 O)AFFE_MATERIAU materia 
@ MAILLAGE mesh 
SO are 
@ TOUT oul 
© MATER : material | 
4_| Expanations: After starting the command file with the DEBUT command, we defined the material giving the 
material characteristics (Young’s modulus and Poisson’s Ratio). Then we specified the format of the mesh 
(MED). Next command defines where the mesh is going to be applied (in the group of elements named “Load”). 
Next figure starts with AFFE_MODELE, in this command we assign to the model, the kind of phenomenon we 
want to use, which will be mechanical and the type of modelization, POU_D_T (which is used for the 
Timoshenko beam theory), AFFE_MATERIAU assigns the mesh to the whole (keyword TOUT) material. 
# Point of Reference: In the previous Beam System cases, after AFFE_MATERIAU, we used to assign the Cross 
Section to the Beam. 
# In this case, we are asked to use I Cross-Section. 
# Unfortunately, Code_Aster does not have a command that will give us the possibility to automatically adjust the I 
Section to the beam (C_A has Rectangular and Circular Cross Sections only). 
# Code_Aster has Rectangular, Circular Cross Section options. 
# But it also has an option called GENERALE, which gives you the ability to manually create whatever type of 


cross section you need. 
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Description 


Figure 


AFFE_CARA_ELEM is the command we use to assign our Cross 
Section to the Beam. 

SECTION: GENERALE, is the demanded point. 

GENERALE option though, needs some parameters in order to 
assign the correct cross section. 

The parameters that are needed for our Timoshenko beam 
(therefore “POU_D_T”) are the following: 

A> Area. 

TY > Principal Geometrical Moment of Inertia compared to GZ. 
IZ > Principal Geometrical Moment of Inertia compared to GY. 
JX > Torsion Constant. 

RY > Distance from furthest left point of beam section to the 
centroid. 

RZ > Distance from bottom of beam section to the centroid. 
AY > Shear Coefficient of GY. 

AZ > Shear Coefficient of GZ. 


0 ECAR 
‘© ose 


canna.) 
CALWLARURA GAR 


Twill show you now, how to add Parameters in Eficas 
MultiCatalogue, and afterwards we will see how what we haye to 
do, to acquire those values, 

If you right click on the last command you used, 
AFFE_MATERIAU in this case, you will see a pop up window 
like the one shown in the figure on the right. 

Go to Parametre > apres (=after). 

This means that Eficas will add a parameter, right after the 
AFFE_MATERIAU command. 


These are all the values for each one of the parameters that the I 
Cross Section requires. 


m AFFE MATERIA materau 
PARAMETRE Tyreasa333 
PARAMETRE 

1 PARAMETRE 

1 PARAMETRE Aye2.19255 

1 PARAMETRE 70073 

1 PARAMETRE RVSTS 

PARAMETRE 


Point: In order to acquire these values with Code_Aster, we will need to use a Macro command. In the next step 


we will show you how to do that. 
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3! 
Step $: Creating a Macro Command 


# In order to acquire the parameters for the I Cross Section, we will need to create a macro command, as mentioned 
before. 
# The way this macro command works is the following: 
© We will make a sketch of a face for our I Section. 
© We will mesh that face. 
© We will create a command file, which after running the analysis will provide us (in the .resu file) all the 
necessary parameters, 
# Let's see that in action, open another Salome _Meca study case, save the case in another folder, activate Geometry 
‘Module from the drop down list and sketch the I Section, 


# |Description Figure 
Z._ | We need to create the face in XY Plane. : ee anacon = 
coma stm 
Go to New Entity > Basic > 2D Sketch. foes crane ten =I sees | 
(imme 
Leave everything Default, Coordinate system, Element Type, ——— 
Destination, Point. oe 79 ce 


Input the following Values: 


100,12.5) 
100,102.5) 
(0,102.5) 
12=(0,12.5) 


(0,0) 
(6012.5) 

= (100,115) 
P_10=(40,102.5) 


Finish with “Sketch Closure”. 


2._| Creating a Face out of our Sketch. 
Go to New Entity > Build > Face. 


Input Sketch_1 for “Objects”. 


Face creation frm wites and/or edges 


1 Thy to crete a planar face 


Finish with “Apply and Close”. 
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# |Description Figure ‘| 
3. | Next and final step, will be to create our groups, o er crap = 
For this macro, we are going to need only one group, which will |] 
include the lines of the perimeter. ee a ids 
a 
Go to New Entity > Groups > Create. sone [ST 
a shape ana Subshapes 
Select the 3" Algorithm for Shape Type of Lines and click on manshops @|[er 
“Select All”, since we don’t have any other lines, i Shape Slt etn Son any seated 
ee ae 
1 ony Subahapes ofthe seco sha | Show alsvbshops] 
second shape _@ |[ 
 is36 ee 
ici Renae 
Finish with “Apply and Close". | {napivandciose] _a00'y_| res] isto | 
4._| These steps conclude our Geometry, save in a proper location and move to Mesh Module. | 
5. | Open Mesh Module from the drop down list. > Et meshysub-mesh +x 
[mame [@ [ech 2 
Go to Mesh > Create Mesh. on _————— 
Select Face_I for Geometry. =o 0) cal et 
Input Netgen 1D-2D for Algorithm. eee rt ete ADE2D ‘| 
Input NETGEN 2D Parameters for Hypothesis. Hypothesis [NETGEN 2D Paremeters =] al 3 
Max. Size = 10. = 
Min, Size = 9. 5 nergen 30 2d, Hypothesis [=> a] aa | 
eae [aca Assign a set of hypotheses s 
Finish with “Apply and x sie fie || lie (aa | 
Close”. rae =! py Aoply Help 
6. | Right click on Mesh_I from the Object Browser > - Mesh computation succeed + 
Compute. compute mest $$ 
ca 
Your Mesh information should be the same as shown Name 
in the figure on the right. nee 
Mes Toe 
Total Linear uactatic 
odes : * 
OD Elements: 0 
Balle: ° 
Edges : 58 2 ° 
Faces: 7 7 ° 
Tangles 18 7 ° 
Quadrangles: 0 ° ° 
Polygons ° 
Volumes: 0 ° ° 
verrahedrons: 0 ° ° 
Herahedrons: 0 ° ° 
Pyramids: ° ry ° 
Prine: ° ° ° 
Hexagonal prisms: 0 
Pobhedrons: 0 
Click “Close”. = 
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# |Description [Figure ‘| 

7. | The figure on the right shows the commands that macrocomm | 
conclude the macro command file. Enna (concopintair | 

= © macro.comm 
Starting with DEBUT, assigning MED format to the pre 
mesh, defining mesh to the whole geometry. @ FORMAT: MED 
© b format med 
a : mezh 
MACR_CARA_POUTRE > Macro command file. @ waniace: mesh 
This will create a table in the .resu file with our 
parameters. 
: ORIG INER : 
ORIG_INER > Original axis orientation. INFO: 
GROUP_MA_BORD > Define the borders of the Hearne 
eometry with a group. @ bom” 
i . on & ® b-gma bord 
. NOEUD : naz 
In the part “NOEUD” you must assign a node from Lm FIN 
your mesh model, it does not matter which node you 
will chose. 

8. _| Activate Aster Module from the drop down list. a —- D 
Go to Aster > Add Study Case. Fs sn 
Input the following data, and run the analysi : 

(ome) Se) 
Sem (EE Tae [STH a 
ctr T]he ang 

9, | After running the analysis you will see that nothing happens to Salome_Meca Object Browser. Don’t get 
frustrated, go to the folder where the -hdf, comm, .mess files are located, and track down .resu file. 

Open it, and you should be able to see the following results, as shown in the figure below. 

These are all the parameters Code_Aster produced for us regarding the I Cross Section. 

Apparently, you can see by yourself that we do not need to use all of them (depends on our MODELISATION 
type). a 

DO, | #4STER 11.3.0 CONCEPT nace CALCULE LE 08/04/2016 A 13:54:39 DE TYPE TABLE SDASTER 

uy A Wee) Sie ieee i iw) SRE RE Iii Aint.) ie ie: ole 

S00 4A. TATE INKED TOAD SME 3 7OSTIEAD OSHS SSE <TR ANSE 5 SOME 5 NEAT 9 READ 
# This was the procedure we use to acquire the necessary parameters for the I Cross Section. 
# Save this Macro Study and head back to our Beam Element case, 


mca 


_ VII: LEVEL_OF_GEOMETRY 


1D - 2D - 3D_ANALYSIS_COMPARISON 


Description 


Now that we assigned the Parameter Iso assigned the 
GENERALE Cross Section to our Beam, it is time for the Degress 
of Freedom and the Loads. 


AFFE_CHAR_MECA > Assigns characteristics to our model. 


Fixed Support to the group of Nodes named “Fix”. 
Displacement Support to the group of Nodes named “Disp” 
Load of 1000 N/mm to the group of Elements named “Loa 


1.000 N/mm and not 10MPa [given data] ? 


The equation for the Pressure equals with P= F/., 

where F = Force and § = Area of the surface where the Force is 
applied. [F = P*S] 

The equation for the Linear Pressure equals with q = F/, , where 
F = Force and | = Length of the Beam. [F = q*l] 


In this study case we want the same Forces, 
consequently P*§ = q*1=> q=P * s), (0) 
S = hw [h=height, w=width] 


@ > q=10* 60000; 1000 N/mm 


Next commands will be the followings’ 


MECA_STATIQUE > Linear Static Mechanics with 
superposition of loads. 


CALC_CHAMP > Calculating the elements/nodes results. 


IMPR_RESU > Printing the results which we got from the 
CALC_CHAMP command. 


AFAR: 
ea segue 
© chat aren 


8 cucamir 
epee 
7B tna 
borane 
teste 
eran 


(AIYIAYAZ)AR 
Lek yA) ARERZ) 


Dela 


We are now done with the command file. 


Go to Fishier > Enregistrer, to save the configuration we made, and close Eficas MultiCatalogue. 
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#_| Description Figure 
‘Activate Aster Module from the drop down list. = eseapiaion % 
soy can 
Go to Aster > Add Study Case. ee 
ae : Coe we fron at] fe werersere.cnn 
Fill the gaps similar with the figure on your right. 
Locate your command file from your computer disk. em) 
Choose Mesh_l as your Mesh. ES 
Ce ne fort al nena ow op 
Change to STA11 for Aster Version vase 
‘imeney oa) [TT —— re ma 
Modify the Solver parameters. ee ee 
ri ia 
Finish with “Ok” se 
‘Your Object Browser should have the same looks as the figure | Object Browser ax 
on the right. 
Right click on Solids > Run. @ Update mesh 
Your analysis will start running and when it finishes Post-Pro 
Branch should be visible. eS 
+ 72600,00 SP 
The analysis may take a while, since we used small size of 7 Line a ret 
elements and our solid model is a little bit bigger than the Perle 
analysis we did so far. ° -7-600.0.0 Kal Cony 
alysis so far. ee pails . 
Sep a 


Export to ASTK 


= Hypotheses 
i Apr: Peres FS 
eiicas, Expand All 


Find 
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Step 11: Post Processing the Results 


Description 


The steps to acquire the deformation and the stress result visualization are always the same, and you should be 
familiar with them until now. 


Right-click on 0,- > Deformed Shape > Tick Magnitude Coloring and change the Scale Factor to a bigger 
number, so you will be able to see the deformation more clear 

Right-click on Def.Shape > Show Only > Representation > Surface 

And you should have something similar to the figure below 


Min. Deformation = 3.6e-21 mm, Max. Deformation = 4.65 mm 


For the stress results: 
Result_SIPO_NOEU > Scalar Map > Mode: SMFZ (Stress due to bending Moment MZ) > Show Only. 


We can see that the direct stress has a range from -3.84 MPa to 767.9 MPa. 
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Solid_Elements 
Step 3: Model Geometry 


# Open up Salome_Meca, activate Geometry Module from the drop down list, save your study to a desired location 
and get ready to start modeling. 


# |Description Figure 

J. | To create our Solid model, we will first create a face of our I Section ‘ iach ceeen = 
and then we are going to Extrude it. cooromaeesymem 
Go to New Entity > Basic > 2D Sketch. (aaa ——— 


Leave everything Default, Coordinate system, Element Type, 9 ca 


Destination, Point. 


Input the following Values: 


P_1=(0,0) P_2=(100,0) P_3 =(100,12.5) |e este 
P_4=(60,12.5) P_5=(60,102.5)  P_6 =(100,102.5) [ese eee 
P_7=(100,115)  P_8=(0,115) P_9 = (0,102.5) 

P_10=(40,102.5) — P_11 = (40,12.5) P_12=(0,12.5) 


“o 2 so aera 
Finish with “Sketch Closure” ||| eos | stachcosne| cases | _ so 


2._| Creating a Face out of our Sketch. 
Go to New Entity > Build > Face. 


Input Sketch_1 for “Objects”. 


Face crestion frm wices and/or ages 


| 
a 


IF ny rmcreate planar face 


eres] env | om | 


Finish with “Apply and Close”. 


3. | Extruding our Face. > ‘onsracon by Foon = 
a 
Go to New Entity > Generation > Extrusion. ro ce 
——————— 
Input the values according to the figure on the right. Ee 
$= 
| 
one, ica 
weight: food —~—s. 
Fan ovectene 
Prete 


scale he face oppest tothe base 


stakes: [ES 


Finish with “Apply and Close”. 


seplyand cose] ny | che a 
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# |Description Figure 
4,_| Final step, will be to create our groups. ~ rate Group x 
hap yp@¢ $$$ 
Go to New Entity > Groups > Create. re . ea ce 
‘croup Nae — 
Select the 3" Algorithm for Shape Type of Faces and “Add” | 


Fixed Support, Displacement Support, Loading and Guide 
group for the whole Solid Model (Select All). 


SS 
man shpe @[[Emust 
| Main shape Selection restriction See naaaeeal 
Ne vesteton 
© Geometrical pans of the Second shape 
© only Sub-shapes ofthe second shape 


i 


Hide selected 


Show al sub shapes 


Finish with “Apply and Close”. || | xstzedtse] anos | se _| __ toe 


5._| These steps conclude our Geometry, save in a proper location and move to Mesh Module, 


cumcnar 
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Click “Close”. 


#_| Description Figure 
7. | Open Mesh Module from the drop down list. Ec metvsabmanh 7 
nome [@ [we 
Go to Mesh > Create Mesh. Geometyy  [Extuson? 
Select Extrusion_| for Geometry. 30 | 20 | 20 | 00] 
fgorthm —_[feaen 1030 =] 
Input Netgen 1D-2D-3D for Algorithm, ypothess —[RETCENSD Pomc =] aif af 
Input NETGEN 3D Parameters for Hypothesis. ot $10 
Max. Size = 10. 
Min. Size = 9. % wetgen30 assign a set of hypotheses | 
arguments | Local ss | avply.and Cosel apoly | ___ cose |] elo 
Finish with “Apply and Close”. yan siee | 
tn, ive | 
2,_| Right click on Mesh_I from the Object Browser > Compute. |]~ tuesh computation succeed al 
Compute mes 
Right click again on Mesh_1 > Convert to/from quadratic, to os 
change our elements from Linear to Quadratic form. 
Right click again on Mesh_1 > Compute. 
tear undratie 
Your Mesh information should be the same as shown in the 
figure on the right. 
° ° 
° ine 
° ° 
° ° 


yyy Cuaprer_VIL: Le 


L_OF_GEOMETRY 


1D - 2D - 3D_ANALYSIS_COMPARISON 


Step 5, 6, 7, 8, 9, 10: Creating the command file, run the analysis 


# Activate Eficas Module from the drop down list and start programming. 
# |Description Figure 
1. | Goto Eficas > Eficas MultiCatalogue. 
In the pop-up window, click on Create New or Nouveau. 
Choose ASTER > “Ok”, STAI1 > “Ok” and you have set up Eficas and it is ready for modification. 
Before starting editing, save the command file in the same location where the geometry/mesb/.hdf files are 
currently located, with an appropriate name. 
Fishier > Enregistrer sous. 
In the figure below, you can see the completed command file, explanations will follow. 
2. 1 @ solids.comm (> @ AFFE_CHAR_MECA: Displac a 
Deaur: F@ MODELE! model 
EPL MATERIAU : material & > BoLimeo 
eu: ©. © GROUP NO: Disp 
ee: 200000 be o 
eu os @ m AFFE chan MECA toad 
URE MAILLAGE fresh @ mopeLe model 
S rorwar: neo 
$ btormae mea Load 
Soot maictace mesh 10 
Sonaltace mesh Ou result 
5 @ OFIE PEaU 30 © mooie model 
@'cRour-Ma: Load materiau 
@e a arre moveve model 
@ wantace mest 
8 @ AED race 
© Tour our 
PHENOMENE BECANIQUE Dieplac 
MODELISAnON =o Load 
©e wm aree wareRiau Materiay — @) Gm caLc came result 
6 waliace mesh @ nesutrar result 
8} Ares $e 
S tare ‘material ae 
Oe a acre can meca: Fed Oem SiGM Noe 
@ MODELE model @ FORMAT. MED 
& @ DOLIMPO & bformat, med 
Fx 3 threvent 
6 aS Rese 
8 © MALLAGE mesh 
8 © kesuttar result 
Mamie coe Bisplac i 
 nobeLe model 
Disp 6 
ac Se homcnan (OEPL StGM_NOEU) 
(©) @ PRES _REP aia Py bee 
‘GROUP_MA Load s aN? ee 
Ce pres Fr} 
3._| Starting with DEBUT command. 6, AFFE_CHAR_MECA > Applying mechanical 
1, DEFI_MATERIAU > Defining the mate Supports, Fixed Support. 
Structural Steel in this case. 7, AFFE_CHAR_MECA > Applying mechanical 
2. LIRE_MAILLAGE > Reading the mesh, MED Supports, Displacement Support. 
format. 8, AFFE_CHAR_MECA > Applying mechanical 
3. MODIL_MAILLAGE > Modifying the mesh, Loads, Pressure of 10 MPa. 
orientation of the normals of the skin of 3D elements. _| 9. MECA_STATIQUE > Linear Static Mechanics 
4, AFFE_MODELE > Applying characteristics to the | with superposition of supports and loads. 
model, Mechanical phenomenon, and 3D modelization, | 10. CALC_CHAMP > Calculating the elements/nodes 
5. AFFE_MATERIAU > Applying the characteristics | results. 
of the material to the mesh model. 11, IMPR_RESU > Writing down, printing the 
results. 
Finish with FIN Command. 
# You will notice that in the section where we applied the supports, we didn’t use rotation restrictions. That is 


because when it comes to solid elements there is no rotation, In Salome_Meca when you apply supports or loads, 
the modifications are applied to each node separately, so when you Fix a node in the top left part of the geometry 
and one more node in the bottom left part with DX, DY, DZ = 0 the model would not be able to rotate, That is 


why Rotation restrictions are redundant here. 
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Description 


Figure 


We are now done with the command file. 


Go to Fishier > Enregistrer, to save the configuration we made, and close Eficas MultiCatalogue. 


Activate Aster Module from the drop down list. 

Go to Aster > Add Study Case. 

Fill the gaps similar with the figure on your right. 
Locate your command file from your computer disk, 
Choose Mesh_I as your Mesh. 

Change to STAI for Aster Version 


Modify the Solver prameters. 
Finish with “Ok” 


pene B55 
Gemmmane [roves 2] af fra peaen Mosca 


ey [Fomcuec owes] gt 


sever [eorot —s] aterverion [ST 5] pred seers 
OS oe 


eo =f 
Stee F 


a 


Your Object Browser should have the same looks as 
the figure on the right. 


Right click on Solids > Run. 


Your analysis will start running and when it finishes 
Post-Pro Branch should be visible. 


The analysis may take a while, since we used small 
size of elements and our solid model is a little bit 
bigger than the analysis we did so far. 


mah 


viyptheses 
£ iitios' & Expat toasre 


Fd ‘Gee 
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ip 11: Post Processing the Results 


# Activate Post-Pro Module and start interrogating the results of our Finite Element Analy: 


# 


Description Figure 


The steps to acquire the deformation and the stress result visualization are always the same, and you should be 
familiar with them until now 


Right-click on 0,- > Deformed Shape > Tick Magnitude Coloring and change the Scale Factor to a bigger 
number, so you will be able to see the deformation more clear. 

Right-click on Def.Shape > Show Only > Representation > Surfaceframe. 

And you should have something similar to the figure below. 


Min. Deformation = 4.2e-20 mm, Max. Deform: 4.79 mm 


‘Same steps for the Normal Stress 
767.602 


Right-click on 0,- > Scalar Map > Change Scalar Mode to SIXX. 
Right-click on Scalar Map > Show Only > Representation > Surface. 
And you should have something similar to the figure on the right. 309.917 


Min. Normal Stres 
Max. Normal Stré 


43.1 MPa 
767.6 MPa 12.2315 


-365.454 


-743.139 
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iii. Surface_Elements 


Step 3: Model Geometry 


# Open up Salome_Meca, activate Geometry Module from the drop down list, save your study to a desired location 
and get ready to start modeling. 
# Unfortunately Salome_Meca does not have an automated option to acquire the mid-surfaces out of our solid 
model, like ANSYS Workbench or other software. We will have to sketch down the Surface model by ourselves. 
# Description Figure 
T._| In order to create our Shell model, first we will need to create the ~ Feit conacton ry 
middle lines of every face, and then Extrude the lines. ———— 
Go to New Entity > Basic > Point. 
Input the following points: 
P_1 =(0,6.25) P_2=(50,6.25) P_3 =(100,6.25) 
P_4= (0,108.75) P_5 = (50,108.75) P_6 = (100,108.75) 
Finish with “Apply and Close”. soviy_| cose | __ Ho 
2._| Next step is to connect the Points with Lines, 
Be careful here with the points because it can be a little complicated. ca 
Go to New Entity > Basic > Line. 
Connect the following Points: aa 
P_1 with P_3. Click “Apply” int ¢ | [= | 
P_4with P_6. Click “Apply” vance [SO 
Se 2 fe [acess 
P_2 with P_S. = 
Finish with “Apply and Close”. ‘iwrdci] won | come | vee | 
4_| Extruding the Lines. = Conscious * 
ern tS 
Go to New Entity > Generation > Extrusion. [a 
—__  __—— 
| 
Fill in the gaps according to the figure on the right. [Rese EEE 
And do the same for Line_2 and Line_3. Be Shapes + cor = 
CC | 
vee [e]foe 
| 
| Fath tne 
neve 
TF Sealethe ce ppt th boe 
——— Or 
L 
Finish with “Apply and = == 
piri | Aectncw ton | 900 | wie |] 
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# | Description Figure 1 
5. | After creating the 3 different extrusions from the lines, we 
need to somehow connect them. Fuse command will help us | | + 
with that, oy 
: _—— 
Go to Operation > Boolean > Fuse. OO 
Here is the tricky part, Fuse commands gives allow us to use || -#94mnss 
only 2 arguments to Fuse, and we need 3 (3 lines, 3 object _@ | [Ena 
extrusions). To make it work, first we will chose 2 out of 3 Sen ella 
extrusions, like the figure on the right. 
Finish with “Apply and Close” coeisence 
Go to Operation > Boolean > Fuse again, but this time for 17 Set presentation parameters and sub-shaps frm arguments 
Object | input the “Fuse_1” we created before, and for I sd prefix to names of restored sub-shapes 
Object_2 input the extrusion you didn’t use before, 
Finish with “Apply and Close” sestyand.cose| app | _cise_|_ ver 
6. | Last step is creating our groups. 
Shape ype $$ 
Go to New Entity > Group > Create. = « a ce 
Graup ame 
Select the 2" Algorithm for Shape Type Lines. — — | 
Fixed Support: Add the lines where we want our fixed support | [Hoc ggg 
to be located. You will notice that there are 5 lines instead of 3, si x 
that would not be a problem. Hen share @ | [re 2 
Displacement Support: Add the lines where we want our foe Shon only sect | 
displacement support to be located, opposite side of fixed ileee! Hide wee 
ene © Geometrcl pas ofthe Second shape | ESS 
PPOM © Only Subshapesofthe Second Shape Show allsubshapes] 
Select the 3" Algorithm for Shape Type Faces. esceer (| 
Loads: Add the top faces where we need to apply our pressure. sect 
Guide: Select All, to add everything in the Guide group. ai 
serve 
‘ep we Ce 
7. | In this study case, we are going to need two more groups, i= ‘dit Group %)| 
which will be describing the different thicknesses we got. Shope ype 
c- ¢ cau ce 
Selecting the 3% Algorithm again for Shape Type Faces. I ——E——— 
Group Name: TpBtThie (Top and Bottom Thickness) = 
Add the top and bottom faces which have thickness of 
Hain Snape Ar Sub shapes 
12.5mm. 
ee | 
Group Name: MidThic (Middle Thickness) pores esau Sieve 
Add the middle face which has thickness of 2 Biel Hide saacted 
© Geometrical pats ofthe Second Shape | ———— 
© Only Subsshapes of the Second Shape | Snew Mi sub-shapes| 
second shape _@ | 
al Select All 
Fa mH 
bosiy coe _| Hep 
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Step 4: Meshing Geometry 
#_[ Description Figure 
J. | Open Mesh Module from the drop down list. > Create mesh aa | 


Go to Mesh > Create Mesh. 


Select Fuse_2 for Geometry. 


falgorthm ——_fretgen 1b25 >] 
Input Netgen 1D-2D-3D for Algorithm. typothess [RETENZO anes =] al Bl 
Input NETGEN 3D Parameters for Hypothesis. —— Tl 
Max. Size = 10. 
Min. Size = 9. 4 wergen 30 assign a set af hypotheses 4 
arguments | toca ss | ‘nly and cess] sopy | chose | _—Hee 
Finish with “Apply and Close Max. Size | 
Min. sie | 


Name [Wea 
‘Geometry [@ [Fuse 2 


30] 1} oo] 


a ‘Mesh computation succeed 
‘Compute mech 


©@ 


Name: 


Mech 2 


Fiangles 
Quadrengles : 
Polygons, 


Volumes : 
Tetrahedrons : 
Hexahedrons : 
Pyramids : 

Prins 
Hexagonal prisms 
Polyhedrans « 


cummcear 
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Step 5, 6, 7, 8, 9, 10: Creating the command file, run the analysis 


# 


Activate Eficas Module from the drop down list and start programming. 


# 


Description Figure 


Go to Eficas > Eficas MultiCatalogue. 

In the pop-up window, click on Create New or Nouveau. 

Choose ASTER “Ok”, STAI1 > “Ok” and you have set up Eficas and it is ready for modification. 
Before starting editing, save the command file in the same location where the geometry/mesby/.hdf fil 
currently located, with an appropriate name. 

Fishier > Enregistrer sous. 

In the figure below, you can see the completed command file, explanations will follow. 


sracecomen yA ATFE CHARUECA: Supas BL IMPRAESU 


BUT. 8 moose rode! Fon ep 
IRE MALAGE ¥ © voto: @ byfermat med 
FoRwat  @ DL mn © besten 
4 berm. med © Grou Fe + nis 
+m ber. cnour mesh © ox: a ¥ @ na: 
© waiace mech, © bv. o © aLLace mesh 
© 4 CHEA GHOUP MA. © or ° (© esutr rest 
@ wow ou + @ ooLiro.2 $ binfo.mes 
your out © GROUP No: Dip 
¥ © cata crour no: 8 oz ° 
¥ @ boption BL AFFE CHARMECA: toade 
(© TOUT. SROUF_MA out (8 Moose rogel (GEFLSiEQ NOEL) 
1B MooL MAILAGE mesh ¥ © pacer. 
8 MALLAGE rest © scnour sa oad 
¥ © one Non coque © pres. 0 
© GROUP MA tose sca StATQUE: result 


Fe MODELE rode! © mooeie ‘model mesh 


MARLASE: rest (© chastaaren aera ASU resut 
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Up until now you must have acquired enough knowledge to be able to understand the purpose of each command, 
For that reason we will discuss only the newly added commands. 

# ORIE_NORM_COQUE (MODI_MAILAGE) > Orientation of the normals for the Shell Element. 

# COQUE_1/_2 (AFFE_CARA_ELEM) > Applying thickness to the given groups (we got two difference 
thicknesses in our geometry). 

# ANGLE_REP (AFFE_CARA_ELEM) > Nautic Angles, used to define a local coordinate in the tangent plane 
of a 2D element. 

# POST_CHAMP > In order to visualize Normal Stress results 
(moy - bottom, sup - top), we have to include this command, be 
Equivalent Stress. 


‘SIGM_NOEU) in particular section of thickness 
se DKT Modelization only allows us to prin 


In the Surface elements studies, like in Solid elements. That means it is not nec 
because, by applying to all the elements, of the groups assigned, DX. DY. DZ 
rotations to be equal to 0. 


y to add rotational restrictions 
, you automatically force the 
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1D - 2D - 3D_ANALYSIS_COMPARISON 


Description 


Figure 


We are now done with the command file. 


Go to Fishier > Enregistrer, to save the configuration we made, and close Eficas MultiCatalogue. 


‘Activate Aster Module from the drop down list. 

Go to Aster > Add Study Case, 

Fill the gaps similar with the figure on your right. 
Locate your command file from your computer disk. 
Choose Mesh_lI as your Mesh. 

Change to STAI for Aster Version 


Modify the Solver prameters. 
Finish with “Ok” 


‘gesuboplistian 7 


study cas defo 


‘command ie [Fam dak =] _ | foracesurcesurice corm 


es [Frc a __e | [ren 


[fecarost x] ate verion [ST — 3] afc serves 


fxeeuton mase TEE] interactive ohow up 


Sele pramters 


‘etal memory ned [7028 a] Time | 
counumber FT ____] Save rut database FP 


‘Your Object Browser should have the same looks as the figure on the right. 


Right click on Solids > Run. 


‘Your analysis will start running and when it finishes Post-Pro Branch should be 


visible, 


The analysis may take a while, since we used small size of elements and our solid 
model is a little bit bigger than the analysis we did so far. 


2 Spt toa 
spin at 
Fo cue 
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1D - 2D - 3D_ANALYSIS_COMPARISON 


Step 11: Post Processing the Results 


Activate Post-Pro Module and start interrogating the results of our Finite Element Analy: 


# 
# |Description Figure 
1._| Right-click on 0 > Deformed Shape > Tick Magnitude Coloring and change the Scale Factor to a bigger 
number, so you will be able to see the deformation more clear. 
Right-click on Def.Shape > Show Only > Representation > Surfaceframe. 
And you should have something similar to the figure below. 
Min, Deformation = 6.98e-21 mm, Max, Deformation = 4.76 mm 
2. 
3. | Same steps for the Normal Stress results. 
Right-click on 0,- > Scalar Map > Change Scalar Mode to SIXX. 
Right-click on Scalar Map > Show Only > Representation > Surface. 
And you should have something similar to the figure on the right. 
Min. Normal Stres: 03 MPa 
Max. Normal Stress = 784.7 MPa 
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Chapter_VII: Tunning Fork 


In this next chapter we are going to get familiar with Modal 
I 


We will study a Tuning Fork and try to acquire the most logical frequenc 


output 
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Step 1: Purpose of the FE Analysis/ Description of the Problem 


In this chapter, we want to perform a Modal Analysis to investigate the natural frequencies of a Tuning Fork. The 
specific tuning fork is designed to tune chamber A 440Hz. In the case that the tuner does not meet the stated requirements, 
we will modify the geometry, material or mass of the tuner in order to get the correct frequency output. 

For that reason, we will set from the start some parameters (Parametric Model) which will help us to modify the 
geometry’s dimensions easier and without the need to sketch the tuning fork from scratch. 


Step 2: Input Values for the FE Analysis 


Initial Dimensions Material 
a= 2.5mm; d=8mm; Model Material: Homogeneous, isotropic and linear elastic continuum. 
20 mm; 1 = 100 mm; Structural Steel: Young’s Modulus = 200 GPa; 
5mm; 2mm; Poisson's Ratio =0.3 
7 mm; Copper Alloy: Young’s Modulus = 110 GPa; 


Poisson's Ratio = 0.34 
Supports 


Fixed Support to “b-dimension” line (drawing will help) 
Symmetric model. 


cumcar 
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Step 3: Model Geometry 


# Open up Salome_Meca and right after the GUI loads, activate Geometry Module from the drop down list. 
# On this chapter we want our model to be parametric. 
# There are a couple of ways to create a parametric model in Salome_Meca Platform, but because most of them are 
really complicated, they include coding inside, here we are going to use the most simplest one. 
# | Description Figure 
1. | Having the Geometry Module activated, go to File > | File 
Notebook. D New ma 
Notebook is a very useful application, which allows us to input | gq Open... ctri+0 
as many parameters as we desire. Reopen 
Notebook saves those parameters and we are able to use them | |i) CONRERES: Ge 
in any Module we want. (Geometry, Mesh) X Close Ctrl+w 
Let’s see how it works. tel Save crits 
Save As... Ctrl+Shift+s. 
Dump Study... Ctri+D 
Notebook... Ctri+K 
Load Script. Ctrl+T 
Properties. Ctri+P 
Import. Ctrl+t 
Export. Ctri+eE 
Preferences... Ctrl+R 
Most Recently Used » 
Exit ctri+Q 
2. | On the left column, as obvious, we arrange the Parameters Salome NoteBook +x 
Name, and on the right column we input the Parameter’s rs variable Valve 
Value. 1 
For our geometry I needed the following parameters: z 
a 
o 
b 
c a Update Study 
Finish with “Apply” penton eal Jae Dae 
Leave Notebook open at all times, because we are going to add 
more parameters as we move forward with this study case. 
3. | Salome_Meca Geometry Module, mostly uses coordinates to sketch/ draw the models. Using Notebook for 
sketching is a little bit trickie because we are going to need lots of different parameters 
To avoid fulling up the study with a wide range of parameters, Notebook gives us the options to use Equations to 
express some of the parameters. 


yy Cuarrer_VIL : TUNNING FORK FREQUENCY, PARAMETRIC_MODEL 
# |Description Figure 1 
4, | We have set our parameters. Have in mind that as we move a 20 herch constuction x 
forward this study, we might need to add more parameters in cootlnate system 
order to help us with the geometry/mesh. | [oo 
Go to New Entity > Basic > 2D Sketch. (OEE. Gi 
Element pe 
As you already know for the construction of a 2D sketch, we ee 4 ce 
were pointing out each point’s coordinates (x, y). 
Destination 
Now instead of arranging numbers for Values, we will input (ia 
the parameters we set before. gids Cosy 
paint ‘akon al Pramerers 
For example: The first point, which originally would have & absolute none Hangentan 
been (0, 0) now according to our parameters will be (0, 0). 7 Relative © Radius 
© section ec 
For this geometry we will need 3 more parameters that we 
haven't set it up before, f= a +d (13) ; values 
j=h+e (120) <p - eG 
g= 2 ee 
g=b+e (25) =| 
Next points: (0,0) (0,a) (e,a)_(e, f) 
GH Go we) (go) pan || eae ort sen 
Finish with “Sketch Closure” | 
5. | In order to create the whole model, first we need to create a > create Face x 
face of the 2D sketch. —o—$ i 
ce 
Go to New Entity > Build > Face. 
a 
Finish with “Apply and Close” a = 
[Face creation fram wires andjor edges 
objects @ | [a 


Wy tocreate « planar face. 


toy | gue |_ vow | 


6. | Goto New Entity > Generation > Extrusion, e onstruction by Extrasion 
- eatrason 
Fill in the gaps according to the figure on the right. ory or) 2 
ect name 
Height: c (c=8mm). Ls 


face shapes + Vector 


wwe @ [fe 


Finish with “Apply and Close” 


vector e] f= 
weights fq 
Poth oirections 

T Reverse 

TF scale the face oppeste to the base 
TT | 
Advanced options 

Proview 


sosivand cose] Aoply | cose | ten 
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# 


Description 


‘Now we need to create our fillets, 
Go to Operations > Fillet 3D. 


Choose the 2" Algorithm which will allow us to select edges, 
and input the same values as shown in the figure on the right. 


Finish with “Apply and Close” 


Do the same for the 7mm fillet, but be careful to choose the 
first fillet as “Main Object”. 


Finish with “Apply and Close” 


Figure "| 
ice ce cw 
ss 
rene [Fe 
‘Filet On Selected Edges ——— — 
| vin cojece [[@] [eamaset 

| selected edges [@ [Ensen iay 
. adiea: [Ra] 

je ee | | 
———— 
even 


Now itis time to create the mirror image. 
Go to Operations > Transformation > Mirror Image. 
Input the following values as shown in the figure on the right. 


Finish with “Apply and Close” 


Mirror image was created a little bit away from the location we 
wanted it to be. 


Go to Operations > Transformation > Translation. 
We want our mirror image to move 8mm to the OZ direction. 
Make sure you uncheck the "Create a Copy” 


Finish with “Apply and Close’ 


episniexe|| aeey 


ee oe 


‘cguments = 
| opjcts [@ [roe 
jf £f 
| 
Cs | 


1 i 


as 
1 Set preottn primar and abanaes fom guns 
fete eee seer 
1 Fw 


only and Cose| __Aonly cose | Heb 
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9. 


Last step before arranging our groups, will be to couple the 
extrusion with the mirror image. 


Go to Operations > Boolean > Fuse. 


‘And input the following values as shown in the figure on the 
right. 


Finish with “Apply and Close” 
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# |Description Figure a 
70, | We are going to need only two different groups for this study | ~ fit cro zs 
aaa re 
© 6 cu ca 
Our Fixed Support which would be in the left top and bottom — | | -ereyp ame $$ 
faces. ame [Fx 


‘Main Shape And Sa-hape2 $$ 
Wain shape _@ | [EE 
ain Shape Selection restieton: 
6 Norestcton 
© Geometrical pars ofthe Second Shape 
© Only Subshapes ofthe Second Shape 


Scan shage_@ | 


‘Show only sect 
Hide selected 


‘show all sub-snapes 


select 
And a Guide group which will contain all the faces we created. | | |** 


add 


Finish with “Apply and Close” — 


(Gesirsne ise] _ soe | see] _ veo 


T1__| This concludes our Geometry creation session. Save in a proper location with a proper name. 


Step 4: Meshing Geometry 


# After finishing from the Geometry creation, go to the drop down list and activate Mesh Module, 


# |Description Figure 
7. | Meshing our model, > crentemesh ae 
name [et 
Go to Mesh > Create Mesh. << 
Input for Algorithm: NETGEN 1D-2D-3D 30 | 20 | 1 | op | 
Input for Hypothesis: NETGEN 3D Parameters, with Aeoitim Rewer t2020 
maximum and minimum size of elements 3 and 2. typothess  [RETGEN SO ramees =] | 23) 
Add, Hypothesis [=None> a 
Finish with “Apply and Close” a 
‘assign a set of hypotheses 5 


eae | er 
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2. | Right-click on Mesh_I > Compute. x Mesh computation suceed +x 
Right-click on Mesh_1 > Convert to/from quadratic, 
Right-click on Mesh_1 > Compute. 
Mesh information shown on the right figure. near Ousarate 
° ae 
° ° 
° 200s 
erahedrs: 0 ° ° 
ryamide; ° ° 
rene ° ° ° 
ecagenl prams 0 
Povhedions: 0 
cose 


Step 5, 6, 7, 8, 9, 10: Salome_Meca Modal Analysis Wizard 


# After completing our geometry, arranging the correct groups, using the most suitable algorithm and hypothesis to 
express our mesh, now its time to make some more modification for our modal analysis. 


#_| Description Figure 
1, | Select Aster Module from the drop down list. eS ‘Ge-subapplication x 


Modal analysis 


Go to Aster > Wizards > Modal Analysis. a 
at nd of odd you art wen? 

First, define the type of model. Click “Next”. ETT = 

Input Mesh_1 for mesh, Click “Next”. 

In the next window insert Young’s Modulus = 2e+11, 

Density = 7850 (structural steel) and Poisson’s Ratio = 0.. 

Click “Next” 

Fix Group (DX=0, DY=0, DZ=0). Click “Next” 

Set Number of Modes (N) = 8. 

Save to a proper location, 


Finish with “Finish” 
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‘Our Modal analysis will appear in the Aster branch, 
Right-click on it and choose “Run”. 


After Code_Aster runs the modal analysis, and if everything 
are set correctly, Post-Pro branch will appear, containing all 
the information we need. 


Expand Al 
Find cuit 
When Code_Aster finish running the modal analysis, Post-Pro T@ Coomaty 
vill ape: Bp Mesh 
branch will appear. a 
i: © Algorithms 
In the “Fields” area, there are the 8 different modes (which we | @ 8 ® Mesht 
selected before) showing the different frequencies. B& Aster 
A modatanalysis 
- aaalecn Br  Post-Pro 
‘You can easily notice that the frequency outputs are smaller Ea aiialotia feed 
than expected. MAL 
Families 
P Groups 
That is because we used millimeters when we created the Qed. 
geometry and the rest are in basic SI Units. There is a © MODES_DEPL, - 
consistency table of units that we have to take into 4 ied 
consideration every time we perform a complete Finite 0378659, 
Element Analysis, which is shown below. 0510918, - 
0.525043, - 
2.26755, 
All we need to do now is to change our geometry to meters. eles 
3.a5e61. 
3.4929, ~ ll 
tong | rorce mace | ase, | wrasse | 6 | cup | rorce | sume 
i ™ N Pa ko kare? 10 eso7 | m | on | Pa 
misc 
2 mm | oN | mPa |torma | umm? or 0 sso7_[mm| on | mPa 
Maire! mm/sec? 
3 a | wr per | sus | swore 10 zaiz | x | wr | per 
feeet 
4 in | wr pei ibt br 10 386 iw | we | psi 
secrin | sectin: invecc? 
s im |r pal ir wna | 2sonioa | seca | m | me | pai 
infeee! 
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# [Description [Figure 1 
5. | To change our geometry units from milimeters to : isis . 
meters, all we have to do is duplicate our mesh and use || —seae 
a seale factor to adjust the units. on a 
Activate Mesh Module. Arguments. 
; : me [Mesh | Setter 
Go to Modification > Transformation > Scale a | —— 
wTiusfaina: Geico ane kebis ae 
trom Je |xf0 ayo _azf ai 
First click the “Select whole mesh, sub-mesh or group” || | <eaieractor | 
option. 
© Move Elements 
‘Then fill in the gaps and choose the correct options as eae at F 
shown in the figure on the right. Sasa ses 
[ret saa ——— 
Finish with “Apply and Close” Trevew 
apes] Higa || Tee | 
6. | Right click on the new mesh, “Mesh_I_scaled” and compute it. 
You will notice the the mesh information is not like the others. The thing we did here was, to copy our previous 
mesh with the previous Algorithms and Hypothesis. The only difference is that with the help of the Scale 
Transformation of 0,001 we converted the milimeters to meters. 
7, | Activate Aster Module, right click on the modal analysis and choose “Edit”. 
You just need to change the “Mesh_1" to the new scaled mesh, “Mesh_1_scaled” 
Click "Ok", right click again and “Run” the analysis 
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# |Description 


Figure 


7._| Now that we have the correct units and the frequency outputs seem normal, we can add deformed shape to each 
one of the frequencies to see how much will the deformation be and also in what way the tuning fork is going to 


be deformed. 
2, [2a mare 


r+ moda analysed 
5 fide 
MODES _pEt 


B DeiShape:2 
pb nape:s 
# Dei shapes 
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Step $: Modify Model 


i. Changing the Material 


# The main purpose of this Modal Analysis, is to get familiar with this kind of analysis and also try to adjust our 
model in order to give us a frequency output of 440 Hz. 
# In order to do that we are going to make some adjustments to our model and see how it reacts. 
# First we are going to change the material of our tuning fork. 
# To do that we will need to change the Young’s Modulus, the Poison’s Ration and the Density inputs. 
# We are interested in using Copper Alloy instead of Structural Steel. 
# Copper Alloy’s material properties: Young’s Modulus= —I.lel 1 Pa 
Poison’s Ratio = 0.34 
Density 8300 kg/m* 
# |Description [Figure 
7. | Open your “yourname”.comm file with your creat an 
favourite text editor. a ae Renney 
one oo 
Go to DEFI_MATERIAL lines and change the orcas | eat 
inputs according to the new material’s Freda. pest Unc nrsorcuereo, 
properties, eae 
Save your changes and head back to anaes 
Salome_Meca. 
Activate Aster Module, and re-run the Analysis. 
# When the Post-Pro branch appears again, you will notice that the frequency outputs are changed, but instead of 


getting closer to 440Hz, they went smaller. 
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ii, Changing the Dimensions. 


# Since, changing the material did not help us get the required frequency results, we will try now and change some 
dimensions from our tuner and see how it affects the frequency results. Changing all the dimension would not be 
much of a help for us, we will modify the I=length dimension and possibly the thickness of the tuner. 


# This is why we assigned the dimensions into Parameters (Notebook) in the beginning of our task, 
Now we just have to open up Notebook and change the parameter inputs and just re-run the analysis. 


3 


# Let’ see how we are going to achieve that. 


#_| Description Figure 
1. | Go to File > Notebook 7 Salome NoteBook + 


Change the “h” Parameter, which is the length of the fork, sarisbie we 


from 100mm to 90mm. eid 2s 
2 |o ° 
Click on “Update Study” 3 [> 5 
afc 8 
Then the whole UI of Salome_Meca is going to change and in| | >}, ae 5 
the Object Browser there will only be Geometry and Mesh 
branches. i i 
rh a 
a feo 2 


Remove Update Study| 
soy | com | 


2, | Unfortunately we will need to re-compute the mesh, create a 
new modal analysis wizard and run the analysis in order to see 
the new frequency results. 


The steps you need to do now are: 
* Right click on Mesh_1 and choose compute, right click again 
and convert to quadratic, compute again. 

* Right click on Mesh_1_scaled and choose compute, right 


fg weaacanage 


‘ . % @ Post-Pro 
click again and convert to quadratic, compute again. Enea 
* Activate Aster Module. eM 
Families 
* Go to Aster > Wizards > Modal Analysis. Gare 
* Input the correct data and save to a proper location. Fields 
5-MODES_DEPL,- 
* Remember that we will use Copper Alloy as material here. seat 
* Right click on the new modal-analysis and Run it, -naa36, 
337251, 
; Sagar, 
The figure on the right shows my own results, as7772, 
ieae 13, 
eaurar. 
3042 07. 


3053.08. 


GUI states 
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# You can see that the frequency outputs are a little bit bigger than the previous ones. Which means that as the 
length of the fork gets smaller the frequency gets higher. 


# In order to hit 440Hz, we will need to make the fork’s length a little smaller. 
Open Notebook again and change the length parameter from 90mm to 80mm, 


3% 


# After re-computing everything with the length set to 80mm, my results are: 


'3-moves_oer 


Snla8, 
3815.04. 
> Modal nalts 
GU Sates 


# We will give it one more try to 75mm and see if we get closer to 440Hz. Here are my final issn 
results: 
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